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ABSTRACT 

We have observed 37 Infrared Dark Clouds (IRDCs), containing a total of 159 clumps, in high- 
density molecular tracers at 3 mm using the 22-meter ATNF Mopra Telescope located in Australia. 
After determining kinematic distances, we eliminated clumps that are not located in IRDCs and 
clumps with a separation between them of less than one Mopra beam. Our final sample consists 
of 92 IRDC clumps. The most commonly detected molecular lines are (detection rates higher than 
8%): N 2 H+, HNC, HN 13 C, HCO+, H 13 CO+, HCN, C 2 H, HC 3 N, HNCO, and SiO. We investigate 
the behavior of the different molecular tracers and look for chemical variations as a function of an 
evolutionary sequence based on Spitzer IRAC and MIPS emission. We find that the molecular tracers 
behave differently through the evolutionary sequence and some of them can be used to yield useful 
relative age information. The presence of HNC and N 2 H + lines do not depend on the star formation 
activity. On the other hand, HC3N, HNCO, and SiO are predominantly detected in later stages of 
evolution. Optical depth calculations show that in IRDC clumps the N 2 H + line is optically thin, 
the C 2 H line is moderately optically thick, and HNC and HCO + are optically thick. The HCN 
hypcrfine transitions are blended, and, in addition, show self-absorbed line profiles and extended wing 
emission. These factors combined prevent the use of HCN hyperfine transitions for the calculation 
of physical parameters. Total column densities of the different molecules, except C 2 H, increase with 
the evolutionary stage of the clumps. Molecular abundances increase with the evolutionary stage 
for N 2 H+ and HCO+. The N 2 H+/HCO+ and N 2 H+/HNC abudance ratios act as chemical clocks, 
increasing with the evolution of the clumps. 
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1. INTRODUCTION 

Although far less common than low-mass stars, mas- 
sive stars play a key role in the evolution of the ener- 
getics and chemistry of molecular clouds and galaxies. 
However, the formation of high-mass stars is far less 
clear than their low-mass counterparts for several rea- 
sons. Massive stars are rare and evolve quickly. In addi- 
tion, the regions that host their early stages of formation 
are difficult to observe due to high dust extinction, and, 
with a few exceptions, are located at large distances (> 
3 kpc) . For this reason, most of the research on massive 
star formation has been based on observations of regions 
with current star formation (e.g., "hot cores," Hll re- 
gions), which were initially easier to detect in surveys 
and to follow up in detail. In constrast, objects in the 
earlier "prestellar" or "starless" stage have been much 
harder to find and, in consequence, this stage still re- 
mains poorly understood. 

About a decade ago, however, the regions contain- 
ing the earliest stages of massive star formation were 
identified. Galactic plane surveys revealed thousands 
of dark patche s obscuring the brigh t mid-infr ared back- 
ground (ISO. iPerault et all Il99a MSX. lEgan et all 
199llSimon et al\\2006^ Spitzer. iPeretto fc Fullerll2009l 
Kim et al. .' 120101) . These dark silhouettes were associated 
with molecular and dust emission, indicating they con- 
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sist of dense molecular gas. Such objects were called 
Infrared Dark Clouds (IRDCs). The first studies char- 
acterizing them suggested that they were cold (<25 K), 
massive (~10 2 -10 4 M Q ), and dense (>10 5 cm" 3 ) molec- 
ular clouds with high c olumn densities (>10 23 cm" 2 ) 
(jCarev et all 119981 120001 ) . They also correspond to the 
densest cloud s embedded within giant molecular clouds 
(ISimon et alll2006bft . 

More recent studies on IRDC clumps|3 using dust con- 
tinuum emission, show that the y have typical masses o f 
-120 Mq and sizes -0.5 pc (jRathborne et all 120061 ). 
Spectral energy distributions (SEDs) of several IRDC 
clumps reveal dust temperatures that range between 
16 and 52 K, and luminosities that range from —10- 
10 5 L Q (jRathborne et all 120101 ). Temperatures derived 
using NH3 observations are lower t han dust temper- 
atures and range from —10-20 K dPillai et all 120061 : 
iSakai et alj|2008t iDevine et al]|2011t iRagan et al.H2011ft . 
Star formation activity in IRDC clumps can be inferred 
from high temperatures and high lum inosities as well 
as th e presence of UC H 11 regions (iBattersbv et all 
120101) . hot cor es jRathborne et all 12008ft . embedded 
24 // m sources ([C hambers et all molecular out- 

flows (|Beuther fc Sridharadl2007l: iSanhueza et alll2 010'). 
and maser emission (jWang et all 120061: iChambers et all 
2009). The high masses, densities, and column densities 
of the IRDC clumps, as well as the aforementioned sig- 

4 Throughout this paper, we use the term "clump" to refer to 
a dense object within an IRDC with a size of the order ~1 pc 
and a mass ~10 2 -10 3 M©. We use the term "core" to describe a 
compact, dense object within a clump with a size <0.1 pc and a 
mass <50 Mq. 
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natures of star formation in them, indicate that IRDCs 
are currently forming massive stars. In addition, IRDCs 
harbor numerous candidates for the most elusive earli- 
est phase of high -mass star formation, the "prestellar" or 
"starless" phase (IChambers et al . 2009; Rathb orne et al. 
2010tlRygl et alJl2010HVasvunina et all201lHPillai et al. 
2011 ; iDevine et al.l 1201 lh T Due to the characteristics 
of IRDCs and the large variety of evolutionary stages 
that they harbor, it has been suggested that they are 
the natal sites of all massiv e stars and stellar clusters 
(jRathborne et al.ll2006l 120101 ) . 

With the aim of characterizing the diff erent evolution- 
ary st ages of clumps found in IRDCs, IChambers et al.l 
( 2009) proposed an evolutionary sequence in which "qui- 
escent" clumps evolve into "intermediate" , "active" , and 
"red" clumps. This evolutionary scheme is based on 
the Spitzer/IRAC 3-8 /im colors and the presence or 
absence of Spitzer/MIPS 24 /im point-source emission. 
A clump is called "quiescent" if it contains no IR- 
Spitzer emission (it is IR-dark); "intermediate" if it con- 
tains either an enhanced 4.5 /im emission, the so-called 
"green f uzzies" (also known as E xtended Green Objects 
(EGOs); ICvganowski et aJ] |2008) . or a 24 /im source, but 
not both; "active" if it is associated with a green fuzzy 
and an embedded 24 /im source; and "red" if it is associ- 
ated with bright 8 /im emission, which likely corresponds 
to an H II region. Quiescent clumps are the best candi- 
dates to be in the "prestellar" or "starless" phase, and 
the massive quiescent clumps are the places where it is 
most probable that high-mass stars are in their earli- 
est stages of evolution. An additional category is "blue," 
which describes objects with bright 3.6 /im emission that 
are predominantly unextincted stars. 

Although several properties of IRDCs have been deter- 
mined in the last few years, there is still one that remains 
poorly explored, their chemistry. What is the chemistry 
in IRDC clumps? A re th e evolutionary stages defined 
by IChambers et al.l (|2009f) chemically distinguishable? 
Pres ently, only a few studies have focused on this sub- 
ject dSakai et al |l2008Ll2010Ll2012t [Vasvunina et~all l2011b 
IMiettinen et al.N201lHChen et al.lboill) . 

In this paper, we report initial results from a program 
aimed at better understanding the chemical evolution 
of IRDC clumps, carrying out a multi-line survey at 3 
mm. We observed several molecular lines simultaneously, 
which facilitates comparison between different lines by 
eliminating or reducing observational errors arising from 
uncertainties in telescope pointing and calibrations. The 
main goals of this project are to investigate the behavior 
of the different molecular tracers and look for observable 
changes in the che mistry as a function of th e evolutionary 
stages proposed by IChambers et all (|2009D . 

2. OBSERVATIONS 

2.1. The Sample 

We made single-pointing observations of molecular 
lines at 3 mm using the ATNF Mopra 22 m tele- 
scope0 located in Australia. The source list comprises 
37 IRDCs, containing a total of 159 clumps. This is a 
sub-sample of the 38 IRDCs and 190 clumps identified 

5 The Mopra telescope is part of the Australia Telescope Na- 
tional Facility (ATNF) which is funded by the Commonwealth of 
Australia for operation as a National Facility managed by CSIRO. 



by iRathborne et all (|2006l ) using 1.2 mm dust contin- 
uum emission at 11" angular resolution. The threshold 
limit for a molecular line to be considered as a detec- 
tion was defined at the 3cr level, where a is the rms 
noise in antenna temperature (T rms ; see Table [1]). The 
clumps G015.31 MM4, G027.75 MM5, G028.37 MM5, 
and G030.57 MM4 present no molecular line emission 
above the 3cr detection level; consequently, they are not 
included in any analysis throughout this work. We ex- 
cluded 17 clumps that show signs of star formation at IR 
wavelengths (Spitzer) and are located at a different dis- 
tance than the rest of the clumps within an IRDC (see 
Section 13.41 for more details) . These clumps belong to 
a different molecular cloud that is situated in the line 
of sight. Since the Mopra telescope beam size (38") is 
sometimes comparable to the separation between clumps 
within the same IRDC, ^30% of the sources are con- 
taminated by emission from another adjacent clump. To 
ensure that the observed emission is attributable solely 
to a single clump, and not a neighboring clump, each 
source with an angular separation smaller than the Mo- 
pra beam size from its nearest neighbor was excluded. 
Furthermore, blue clumps are not included in any analy- 
sis because they are likely unextincted stars. Therefore, 
this work will focus on th e remaining 92 IRDC c lumps 
unless stated ot h erwis e. IChambers et al.1 ()2009l ) and 
IRathborne et al.l ( |2010l ) classified the initial sample of 
190 IRDC clumps using a combination of Spitzer/IRAC 
colors and the presence or absence of compact 24 /im 
emission. Their classifications differ in t he method of es- 
timati ng the 24 /im flux measurements. Chambers et al. 

( 2009) used aperture photometry, and IRathborne et al. 

(2010) fitted a Gaussian to the radial profile of the emis- 
sion. Due to this difference, in a few cases 24 /im emission 
is detected for a particular source using one method, but 
not the other. In this work, to avoid ambiguities, if a 
clump presents 24 /im emission from either of the two 
methods, a 24 /im compact source will be assigned to 
the clump. 

Table [2] summarizes information for all 159 IRDC 
clumps: names, coordinates, velocities, distances, dust 
temperatures, Spitzer classification, and if the source is 
an IRDC and if it is used in the present work. 

2.2. Observing Parameters 

The observations were carried out during 2008 July- 
August and 2010 September. They were performed 
with the 3 mm Monolithic Microwave Integrated Circuits 
(MMIC) receiver and the Mopra spectrometer (MOPS)0 
in the "broad-band" mode, resulting in a total bandwidth 
of 8.3 GHz split over four overlapping sub-bands of 2.2 
GHz. Each sub-band has 2 x 8096 channels (1 x 8096 per 
each polarization) with a channel resolution of 0.27 MHz. 
The velocity resolution was ^0.90 km s^ 1 per channel. 
Both polarizations were observed and averaged to im- 
prove the signal to noise ratio. The observed frequencies 
range between 85.540 and 93.840 GHz. This frequency 
range allowed us to observe 15 molecular lines simulta- 
neously (see Table QJ. The system temperatures ranged 
between 165 and 265 K. Typical rms noise values (T rms ) 

6 The University of New South Wales Digital Filter Bank used 
for the observations with the Mopra Telescope was provided with 
support from the Australian Research Council. 
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Summary of Observed Molecular Lines 
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Note. — The critical density was calculated as n cr it = A u i/7 u i, where A u i is the Einstein coefficient and 
7 U 1 is the collisional rate. Values of A u i and 7 u i at 20 K (50 K for SO) were obtaine d for most of the molecules 
from the Leiden Atomic and Molec ular Database (LAMDA ) jSchoicr ct al. 2005). Values of A u i and 7 u i at 
25 K for NH 2 D were obtained from lMachin fc Rouefll (1200611 . 
a Critical density adopted from [L~o et all ]2009f) . 



for each molecular line are shown in Table [T] The on- 
source integration time per object was ^3 minutes. The 
angular resolution of Mopra telescope is 38" and its main 
beam efficiency is 0.5 at 90 GHz[J All the observations 
were performed in position-switching mode with the off- 
position shifted from the target source by 1° in Galactic 
latitude away from the Galactic midplane. The telescope 
pointing was checked by observing nearby SiO masers ev- 
ery ~1 hr, and was maintained to be better than 5". The 
uncertainty in the line intensity varies with the frequency 
from 9% at 88 GHz to 25% at 93 GHz (Jonathan Foster, 
private communication). The initial spectral process- 
ing was done using the ATNF Spectral Analysis Package 
(ASAP) softwareQ The observed molecular lines, transi- 
tions, rest frequencies, upper energy levels, critical den- 
sities, sub-bands, and typical T rms are listed in Table [TJ 
The T rms for each source is given in Table El 

3. RESULTS 

3.1. Line Parameters 

For each molecular transition, the line center velocity, 
line width and intensity of the line were determined from 
Gaussian fits. Molecular lines that show hyperfine struc- 
ture (N 2 H + and C 2 H) were fitted using a multi-Gaussian 
function, with a fixed frequency separation between the 
transitions. The Gaussian fit procedure was carried 

7 See more details about Mopra telescope, receivers, and back- 
ends at http://www.narrabri.atnf.csiro.au/mopra/obsinfo.html 

8 http:/ /www. atnf.csiro.au/computing/software/ 



out i n IDL using the MPFITFUN package (|MarkwardtJ 
2009). Occasionally, two velocity components were de- 
tected in a source. We define the main component, i.e., 
the component associated with the IRDC clump, as the 
velocity component detected in the high-density tracers 
(e.g., N 2 H + and/or H 13 CO + ). The secondary velocity 
component typically corresponds to lower-density gas lo- 
cated in line of sight and they were not used in this work. 
When two velocity components have emission from high- 
density tracers, the brightest was used as the main com- 
ponent. The Gaussian fit parameters of the main velocity 
component for all sources where emission was detected 
(155) are summarized in Table |3] No Gaussian fits were 
carried out for two distinct velocity components with sep- 
arations less than 3 km s _1 or for self-absorbed profiles. 

3.2. Spectra and Detection Rates 

Figure Q] presents the detection rates, at the 3cr detec- 
tion level, of the 10 most frequently detected molecular 
species toward the IRDC clumps. For the N 2 H + and 
C 2 H lines, we show the most often detected transitions: 
JF X F = 123 -)• 012 for N 2 H+ and NJF = l§2->0|l 
for C 2 H. Because the HCN transitions are blended, and 
in addition, exhibit self-absorbed profiles and wing emis- 
sion, the multi-Gaussian fit to the hyperfine structure 
was not reliable. The HCN detection rate was ob- 
tained by inspecting the spectrum for each source by 
eye, instead of comparing the intensity of the Gaus- 
sian fit with the T rms . The CH 3 CN, 13 CS, H 13 CN, SO, 
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and NH2D lines were detected in fewer than 8 clumps; other hand, the molecular line shown in Figure [T] least 
therefore, their detection rates are not presented in Fig- often detected was the SiO J = 2 — 5- 1 with 8 (9%) de- 
ure [TJ The molecular lines most often detected were tections. Detection rates for all detected molecular tran- 
HNC J = 1 — 5- and N2H+ JF\F = 123 — > 012 with sitions are given in Table |U Uncertainties presented in 
90 (98%) and 89 (97%) detections, respectively. On the Figure Q] and Table |4] were determined assuming Poisson 

noise. 



Table 2 

Clumps Properties. 
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MM4 


18:35:19.4 


-07:37:17 


115.0 


7.8 


31.0 


Q 


IRDC 


G024.33+00.11 


MM5 


18:35:33.8 


-07:36:42 


117.2 


7.7 


32.0 


i 


IRDC 


G024.33+00.11 


MM6 


18:35:07.7 


-07:34:33 


114.3 


7.7 




l 


Blended 


G024.33+00.11 


MM7 


18:35:09.8 


-07:39:48 


99.6 


6.3 


23.7 


Q 


IRDC 


G024.33+00.11 


MM8 


18:35:23.4 


-07:37:21 


113.8 


7.8 


31.5 


Q 


IRDC 


G024.33+00.11 


MM9 


18:35:26.5 


-07:36:56 


119.2 


7.7 


52.0 


R 


IRDC 


G024.33+00.11 


MM11 


18:35:05.1 


-07:35:58 


113.3 


7.8 


26.0 


Q 


IRDC 


G024.60+00.08 


MMl 


18:35:41.1 


-07:18:30 


53.4 


3.8 




A 


Blended 


G024.60+00.08 


MM2 


18:35:39.3 


-07:18:51 


115.2 


7.7 


23.0 


1 


IRDC 


G024.60+00.08 


MM3 


18:35:40.2 


-07:18:37 


53.8 


3.8 




1 


Blended 


G024.60+00.08 


MM4 


18:35:35.7 


-07:18:09 


53.2 


3.8 




Q 


Blended 


G025. 04-00. 20 


MMl 


18:38:10.2 


-07:02:34 


64.0 


4.3 




A 


Blended 


G025. 04-00. 20 


MM2 


18:38:17.7 


-07:02:51 


63.5 


4.3 


28.0 


1 


IRDC 


G025. 04-00. 20 


MM3 


18:38:10.2 


-07:02:44 


63.7 


4.3 




Q 


Blended 


G025. 04-00. 20 


MM4 


18:38:13.7 


-07:03:12 


63.8 


4.3 


29.5 


i 


IRDC 


G025. 04-00. 20 


MM5 


18:38:12.0 


-07:02:44 


63.9 


4.3 




Q 


Blended 


G027.75+00.16 


MMl 


18:41:19.9 


-04:32:20 


78.3 


5.0 




A 


Blended 


G027.75+00.16 


MM2 


18:41:33.0 


-04:33:44 


51.4 


3.5 


24.0 


Q 


IRDC 


G027.75+00.16 


MM3 


18:41:16.8 


-04:31:55 


78.7 


5.0 




i 


Blended 


G027.75+00.16 


MM5 


18:41:23.6 


-04:30:42 


79.1 d 


4.8 d 




Q 


No Detection 


G027.94-00.47 


MMl 


18:44:03.6 


-04:38:00 


45.4 


3.1 


44.0 


K 


IRDC 


G027.97-00.42 


MMl 


18:43:52.8 


-04:36:13 


44.6 


3.1 




A 


Blended 


G027.97-00.42 


MM2 


18:43:58.0 


-04:34:24 


19.9 


1.5 




R 


Non-IRDC 


G027.97-00.42 


MM3 


18:43:54.9 


-04:36:08 


45.9 


3.2 




Q 


Blended 


G028. 04-00. 46 


MMl 


18:44:08.5 


-04:33:22 


45.8 


3.2 


38.0 


A 


IRDC 


G028. 08+00.07 


MMl 


18:42:20.3 


-04:16:42 


81.4 


5.2 


22.0 


1 


IRDC 


G028. 10-00. 45 


MMl 


18:44:12.9 


-04:29:45 


47.1 


3.2 




Q 


Blended 


G028. 10-00.45 


MM2 


18:44:14.3 


-04:29:48 


46.9 


3.2 




Q 


Blended 


G028. 23-00. 19 


MMl 


18:43:30.7 


-04:13:12 


80.0 


5.1 


22.0 


Q 


IRDC 


G028. 23-00. 19 


MM2 


18:43:29.0 


-04:12:16 


80.8 


5.1 




I 


Blended 


G028. 23-00. 19 


MM3 


18:43:30.0 


-04:12:33 


80.2 


5.1 




Q 


Blended 
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IRDC Clump Coordinates V; sr a D T dust Class. b Comment' 

Name Number a(J2000) <5(J2000) (km s" 1 ) (kpc) (K) 



0028 28-00 34 


MM1 


1 8-44-1 ^ 


-04:17:54 


48.7 


3.3 




R 


Rl pin rl prl 

D1CI1UCU 


C!028 28-00 34 


MM2 


1S-44-21 


-04:17:37 


84.9 


5.6 




R 


is U11_L± V.J— / V_'' 


C!028 28-fin 34 


1V±1V± vJ 


18:44:13.4 


-04-1 R-Qo 

Urt. IO . UU 


49.7 


3.4 




R 


I -1 1 on 1 1 or 1 

D 1 L 1 1 L. Ll 


0028 28-fin 34 


MM4 


18:44:11.4 


-04:17:22 


48.8 


3.4 


34.5 


A 


IRDC 


C!028 37+00 07 

vJUZjO.O i 1 \l U . \ l I 


MM1 


18-42-52 1 


-03-50-45 

UU .U£7 . rlU 


78.1 


5.0 


33.0 


A 


IRDC 


on28 37+00 07 


MM2 


1 8-42-37 fi 

J.O.rtii. O I . U 


-04-02-05 

Urt.Uii . UU 


80.8 


5.1 


39.0 


x 


IRDC 


nn28 37+00 07 

V~UiiO . U ( | \J U . \) 1 


L\± J.V± vJ 


J_ O . rtO . UO . J. 


-04-06-24 


99.8 


6.4 




R 


Noti-TRDC! 

L\ U11_L± \,i—J V_'' 


nri28 37+00 07 

VJUZjO . U 1 1 UU . U I 


MM4 


18-42-50 7 

.LO.rtii. UU. 1 


-04-03-1 5 

Urt.UO . -LU 


79.3 


5.0 


32.0 


A 


IRDC 


OUzo.o ( H~UU.U ( 


i\ /r a /r " 
1V11V10 


10.4Z.Z0.0 


r\A .ni .Qfi 
-U4.Ul.v5U 


70 fid 
( o.O 


A Qd 
4.0 




it 


No Detection 


OUzo.o / H~UU.U ( 


IVllVLo 


1 Q-/IO-/K1 n 

lo.4Z.4y.U 


HA .flO-OQ 

-U4.Uz.Zv5 


on n 
oU.U 


" 1 

0. 1 


oq n 
zo.U 


A 


1H1JO 


Pfios Q7_i_nn n7 
bUzo.o ( H-UU.U i 


1V11V1 ( 


lo.4Z.0O.v5 


n/i -n7- qi 
-U4.U / .ol 


A« fi 
40.0 


Q O 
O.Z 




It 


IN On-llX-L/O 


OUZo.v5 / H~UU.U ( 


IVllvlo 


lo.4z.4y. ( 


c\a .nn. " /i 
-U4.Uy.04 


1 ri7 
1U ( .Z 


7 " 

/ .0 




T> 

rv. 


INOll-UxUO 


bUzo.o < H~UU.U i 


A/TA/T Q 

iviiviy 


1 Q- A O-A fi 7 
lo.4Z.40. / 


c\a -c\a -ns 

-U4.U4.Uo 


70 A 

( y.4 


" n 
O.U 


OA 
Z4.0 




TDnp 
IHUO 


Pfios Q7_i_nn n7 
bUzo.o * H-UU.U i 


AAA/f 1 n 
1V11V1 1U 


10.4Z.04.U 


c\a -no- Qn 

-U4.UZ.oU 


70 1 
/y.i 


" n 
O.U 




A 


"R 1 on A t±A 

Joicnaca 


pnofi Q7_Lnn n7 
bUzo.o ( -(-UU.U i 


A/TA/T1 1 
1V11V1 1 1 


1 fi- /I O-AO 7 
lo.4Z.4Z. / 


C\A -ni - A A 
-U4.U1 .44 


fin 
ou.y 


1 

0. 1 


Qfi n 
oO.U 


A 
i V 


inuo 


rinoQ 07 i nn H7 
OUzo.o / H-UU.U ( 


1V11V11Z 


lo.4o.Uy.y 


C\A .nfi- "O 
-U4.U0.0z 


Q/1 Q 

o4.o 


" 
z.O 


OQ 7 

zo. / 




lit 1JO 


bUzo.o * H-UU.U i 


1V11V1 lo 


10.4Z.41.O 


c\o. "7. ns 
-Uv5.0 1 -Uo 


77 " 
1 ( .0 


/i 
4.y 


Qn a 

v5U.4 


T 
1 


lltUO 


pnofi "Q nn oc; 

OUZ0.0o~UU.Z0 


1V11V1 1 


ifl./i/i.ia n 

10.44. 10. U 


HQ. "O- Q/1 
-Uv5.0y .v54 


fifi fi 
0O.0 


" " 
0.0 






Blended 


rinoQ "Q nn o" 
OUZo.Oo-UU.zO 


IVllvlZ 


io,/|/(,i c 7 
lo.44. 10. ( 


no. r: n. /i 1 
-Uv5.0y.41 


fi " fi 
oO.O 


" A 

0.4 




A 
1\ 


lolcnaca 


Pnoo "Q nn o" 


IVllvlo 


1 Q. A A . 1 fi n 

lo.44. lo.U 


n/i .nn. /i 
-U4.UU.4o 


Qfi A 
O0.4 


" " 
0.0 


00 n 
zz.U 


r\ 


lltlJO 


\~UZO.OO~uU.ZO 


A/TA/TA 


lo.'l'l. lo.U 


-U'l.UU.UO 


OO.O 


O.'i 




I 


Blended 


O028 o^-00 2o 


1V±1V± u 


18:44:17.0 


-04-02-04 

u^t.u^.u^t 


87.0 


5.7 


30.0 


I 


IRDC 


C4028 53-00 25 


MMfi 

L\± J.V± U 


18:44:17.8 


-04-00-05 

Urt.UU . UU 


86.4 


5.5 




O 
O 


M 1 on r 1 or 1 


C028 53-00 25 

v. I u ^ O • • )- ) uUtiy 


MM7 


18:44:23.7 


-0402 00 


88.6 


5.6 


23.0 


IRDC 


C4028 53-00 25 

VJJUiiO. UO UU- iiU 


MM8 


1 8-44-22 


-04-01 -35 

Urt.U J. .OU 


88.6 


5.6 


28.0 


O 


IRDC 


C4028 53-00 25 

VJJUiiO. UO UU- iiU 


MM9 


1 8-44-1 Q 3 


-03-58-05 


87.1 


5.7 


28.0 


I 


IRDC 


C4028 53-00 25 

VJJUiiO. UO UU- iiU 


MM1 n 


18-44-18 5 


-03-58-43 

UU . uo . , -±u 


86.8 


5.5 


30.0 


O 


IRDC 


C028 67+00 1 3 

V. I W — (. ; . u ' | UU • 1 i ) 


MM1 


18-43-03 1 

J_vJ.^vJ.UvJ. J- 


-03:41:41 


84. 3 C 


5.3 C 


40.0 


R 


IRDC 


C!n28 R7+oo i 3 

VJJUZjO . U I | UU . _LU 


MM2 


1 8-43-07 1 

.LO.rtO. Ul .1 


-03:44:01 


79.2 


5.0 


23.0 


I 


IRDC 


on28 R7+oo i 3 

UUiiu . U I | UU . _LU 


LY± J.V± vJ 


18-42-58 2 


-03-48-20 


103.9 


6.9 




A 


Noti-TRDC! 


VJUOU . 1 i UU • UU 


MM1 


18-46-35 7 


-02-31 03 

Ui-vl . UU 


86.8 


5.5 


27.0 


A 


IRDC 


O030 57-flO 23 


MM1 


18-48-00 

1U-T:U-UU.U 


-02-07-20 

yj j-j . vj 1 -iu 


90.7 


5.8 


34.5 


A 


IRDC 


oo30 57-fin 23 


MM2 


18-47-58 7 


-02-1 5-20 

Ui- -L U . ^U 


111.3 


7.3 




R 


Noti-TRDC! 


oo30 57-fin 23 

VJV'JU . U I UU.iiO 


1V±1V± vJ 


18-47-54 5 


-02:11:15 


95.9 


6.3 


23.0 


I 


IRDC 


nnQfi "7 nn oq 
OUv5U.O /-UU.zo 


1V11V14 


lo.4o.Ul.o 


no- 1 O-Q " 
-Uz. Iz.oO 


Qfi od 
00. z 


" od 
O.z 




C\ 

0, 


No Detection 


PnQn n7 nn i /i 

oUv5u.y ( -uu. 14 


1V11V1 1 


lo.4o.zl.O 


ni . a . 07 
-U1.40.Z f 


77 Q 
( t.O 


" n 
O.U 


Q7 n 
/ .U 


A 
A 


lltlJO 


rinQi no nn in 
OUol.UZ-UU. 1U 


1V11V1 1 


i q. /i con n 
lo.4o.Uy.y 


ni ■ /i " ■ 1 7 
-U1.40. 1 ( 


7fi n 

/ o.y 


/i n 

4.y 




TO 

Jo 


Jorignt rSluc 


rinQi n7 i nn H7 

oUv5i.y * h-uu.u * 


1V11V1 1 


1 Q- /I n. Qfi Q 

Io.4y.v50.v5 


nn- /i " - /i " 
-UU.40.40 


o" n 

yo.u 


fi " 
0.0 




A 


lolcnaca 


riAQi n7 i nn H7 

ouoi.y * h-uu.u r 


IVllvlz 


i q. /i n. Qfi n 
lo.4y.v5o.U 


nn. a fi. 1 fi 
-UU.40. 10 


CiA fi 

y4.o 


fi " 
0.0 




r\ 


Blended 


r^nQi d7j_nn n7 
UUol.y * -(-UU.U i 


IVllvlo 


1 fi- A O- io 'i 
lo.4y.v5Z.o 


nn- /i 7-no 
-UU.4 ( .UZ 


Q/1 1 

y4. 1 


fi A 
0.4 




\4 


Blended 


rinQi n7 i nn H7 

oUv5i.y * h-uu.u * 


1V11V14 


1 Q./in.QQ n 
lo.4y.v5o.U 


nn- /i 7- qq 
-UU.4 / .00 


n" " 

yo.o 


fi fi 
0.0 




T 
1 


Blended 


bUol.y ( -hUU.U i 


1V11V10 


1 fi- /I G-01 o 

io.4y.zi.y 


nn- "n- q " 
-UU.0U.v50 


Ofi " 

yo.o 


fi fi 
0.0 


Q" n 
v50.U 


T 
1 


iituo 


OUv5i.y * -(-UU.U i 


A/TA/Tfi 
1V11V10 


1 fi- A O- Q" n 
lo.4y.v50.U 


-UU.40.44 


O/l fi 

y4.o 


fi " 
0.0 




C~\ 

w, 


Blended 


r^nQi Q7j_nn n7 
bUol.y * H-UU.U i 


A/TA/f 7 
1V11V1 i 


1 fi- /I O-Ofi /I 
lo.4y.Z0.4 


nn-/i q. " /i 

-UU.40. 04 


OQ O 

yo.y 


fi /I 
0.4 


07 n 

Z / .U 




TTO TO 
lltUO 


riAQi o7 i nn ri7 

oUv5i.y * h-uu.u r 


A/TA/TC 
IVllvlo 


i fl./in.on i 
lo.4y.zy. 1 


nn. a ■ 1 
-UU.4o. lz 


CiA " 

y4.o 


fi /I 
0.4 


a 7 n 
4 / .U 


A 
A 


TTO TO/"' 
lltlJO 


ririQQ fin nn m 
OUv5o.oy-UU.Ul 


A/TAT 1 
1V11V1 1 


1 C- "O- "Q Q 

lo.Oz.Oo.o 


UU.4Z.O / 


1U0. / 


7 1 

/ . 1 


o" n 
Zo.U 


Ji 


lltlJO 


ouoo.oy-uu.ui 


A/TA/f 
1V11V1Z 


1 fi- "O- /I 

io.oz.4y.y 


nn-Q7- "7 
UU.vS / .0 / 


1 f\A fi 
1U4.0 


7 1 
( . 1 


/ii n 

41 .U 


"R 

n 


TTO TO ^ 
lltUO 


r^nQQ fio nn m 
OUv5o.oy-UU.Ul 


IVllvlo 


1 q- "o- "n fi 

lo.OZ.OU.o 


nn- Qfi- a q 

UU.v50.4v5 


1 nQ 1 
lUo.l 


7 1 
i . 1 


/iQ n 
4v5.U 


p 
n 


lltUO 


rinQQ fin nn m 
OUv5o.0y-UU.Ul 


A/TAT A 
1V11V14 


1 Q- "O- "fi /I 

lo.Oz.oo.4 


nn. a q .no 
UU.4v5.Uo 


1 nfi 
1U0.Z 


7 1 

r . 1 


oq n 
zo.U 


A 
i\ 


lltlJO 


priQQ fio nn m 
ouv5o.oy-uu.ui 


If ATC 
1V11V10 


1 fi- "O-/I 7 Q 
lO.OZ.4 ( .0 


nn-Qfi- /i 7 
UU.oO.4 / 


in"" 
1U0.0 


7 1 
i . 1 


Q7 n 
/ .U 


A 
J 1 


lltUO 


rinQQ fin nn m 
OUv5o.0y-UU.Ul 


A/TA/T1 1 
1V11V1 1 1 


1 Q- "O- "fi O 

lo.Oz.OO.z 


nn. a 1 . a q 
UU.41.4o 


1 ri7 n 
1U ( .U 


7 1 
/ . 1 


o" n 
zO.U 




TTO T0<^' 

lltlJO 


ririQ/l /I Q I nn OA 

OUo4.4o-|-UU.Z4 


1V11V1 1 


1 Q- "Q- 1 c n 

lo.Oo. lo.U 


ni .o".Ozi 
Ul.z0.z4 


"7 n 

/ .y 


Q Q 
O.O 


qq n 
v5o.U 


A 


TTO T0/"~' 
lltlJO 


OUo4.4o-hUU.z4 


IVllVlz 


1 Q- "Q- 1 Q fi 

lo.Oo. lo.O 


ni -O/i ■ a n 
Ul.z4.4U 


"7 A 

( .4 


Q Q 
O.O 




It 


TO 1 vl ,] „ O 

Jolcnciecl 


OUo4.4oH~UU.Z4 


IVllvlo 


1 Q- "Q-On /I 

lo.Oo.zU.4 


ni ■ 00 ■ oq 
Ul. zo.Zo 


"O /t 

0y.4 


q n 

o.y 




A 
A 


Blended 


ririQ/l A Q I nn O/l 

OUo4.4oH-UU.Z4 


A/TAT A 
1V11V14 


1 Q- "Q- 1 n n 
lo.Oo. ly.U 


ni -OA .no 
Ul.z4.Uo 


"7 fi 

/ .0 


Q Q 

v5.o 




A 


Blended 


ririQ/i /i q i nn OA 
OUo4.4oH~UU.Z4 


A/TA/T " 
IVllVlO 


1 Q- "Q- 1 n c 
lo.Oo. ly.o 


ni ■ oq ■ Qn 
Ul.Zv5.v5U 


"Q n 
Oo.U 


Q Q 
O.O 


oq n 
Zo.U 


A 


TTO T0<"~' 
lltlJO 


ririQ/i /i q i nn OA 
OUo4.4oH~UU.Z4 


A/TAT fi 
IVllVlO 


1 Q- "Q- 1 Q fi 

lo.Oo. lo.O 


ni ■ 07- A 
Ul.z ( .4o 


"Q fi 

Oo.O 


Q Q 
O.O 




T 
1 


Blended 


r^m/i /iQ_i_nn oa 

OUo4.4oH~UU.Z4 


A/TA/f 1 
1V11V1 i 


1 fi- "Q- 1 fi 1 

lo.Oo. lo.v5 


ni -07-1 Q 

Ul.Z / .10 


c;q 1 
Oo. 1 


Q fi 
O.O 


00 " 
zy.o 


T 
1 


TTO TO/"' 
lltUO 


ririQ/i /i q i nn OA 
OUv54.4oH~UU.z4 


IVllVlO 


1 Q- "Q- 1 fi /I 

lo.Oo. Io. 4 


ni ■ ofi.on 
Ul.zO.zU 


"7 O 

/ .z 


Q 7 
v5. ( 


/i q n 
4o.U 


A 
A 


TTO T0<^' 
lltlJO 


r^m/i /iQ_i_nn oa 

OUv54.4oH~UU.Z4 


A/TA/TQ 

iviiviy 


1 Q- "Q- 1 fi /I 

lo.Oo. 10.4 


ni -Ofi- 1 a 

Ul .Zo. 14 


"fi 7 
00. / 


Q fi 
O.O 




r~\ 


Blended 


G034. 77-00. 55 


MM1 


18:56:48.2 


01:18:47 


44.2 


2.9 


40.4 


H 


IRDC 


G034.77-00.55 


MM2 


18:56:50.3 


01:23:16 


42.0 


2.8 




Q 


Blended 


G034.77-00.55 


MM3 


18:56:44.7 


01:20:42 


43.1 


2.8 


30.4 


I 


IRDC 


G034.77-00.55 


MM4 


18:56:48.9 


01:23:34 


42.1 


2.8 




Q 


Blended 


G035.39-00.33 


MM1 


18:56:41.2 


02:09:52 


64.3 


4.2 




R 


Non-IRDC 


G035.39-00.33 


MM2 


18:56:59.2 


02:04:53 


54.0 


3.5 




A 


Non-IRDC 


G035.39-00.33 


MM3 


18:57:05.3 


02:06:29 


55.3 


3.6 




R 


Non-IRDC 


G035.39-00.33 


MM4 


18:57:06.7 


02:08:23 


45.3 


3.0 




A 


Blended 


G035.39-00.33 


MM5 


18:57:08.8 


02:08:09 


90.1 


6.9 




Q 


Blended 


G035.39-00.33 


MM6 


18:57:08.4 


02:09:09 


45.5 


3.0 




A 


Blended 


G035.39-00.33 


MM7 


18:57:08.1 


02:10:50 


45.9 


3.0 


34.0 


A 


IRDC 


G035.59-00.24 


MM1 


18:57:02.3 


02:17:04 


50.4 


3.3 


36.0 


R 


IRDC 


G035.59-00.24 


MM 2 


18:57:07.4 


02:16:14 


45.3 


3.0 


24.0 


A 


IRDC 


G035.59-00.24 


MM3 


18:57:11.6 


02:16:08 


45.0 


3.0 


30.5 


A 


IRDC 



() 
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Table 2 — Continued 



IRDC 


Clump 


Coordinates 


Vi a 

v Isr 


D 


A dust 


Class. b 


Coininciit c 


1 > dUlL 


IN1 1 1 mnor 


cy(J2000) 


<5(J2000) 


( km s ) 










G036.67-00.ll 


MM1 


18:58:39.6 


03:16:16 


53.5 


3.5 


20.0 


Q 


IRDC 


G036.67-00.ll 


MM2 


18:58:35.6 


03:15:06 


53.3 


3.5 


16.0 


Q 


IRDC 


G038. 95-00.47 


MM1 


19:04:07.4 


05:08:48 


42.4 


2.8 


16.0 


I 


IRDC 


G038. 95-00.47 


MM2 


19:04:03.4 


05:07:56 


42.1 


2.8 


40.0 


R 


IRDC 


G038. 95-00.47 


MM3 


19:04:07.4 


05:09:44 


42.4 


2.8 


28.0 


1 


IRDC 


G038. 95-00.47 


MM4 


19:04:00.6 


05:09:06 


42.1 


2.8 


43.0 


A 


IRDC 


G048. 65-00.29 


MM2 


19:21:47.6 


13:49:22 


33.8 


2.5 




Q 


Blended 


G053. 11+00.05 


MM1 


19:29:17.2 


17:56:21 


21.8 


1.8 


38.0 


R 


IRDC 


G053. 11+00.05 


MM2 


19:29:20.2 


17:57:06 


22.9 


1.9 


36.0 


A 


IRDC 


G053. 11+00.05 


MM3 


19:29:00.6 


17:55:11 


2.6 


0.5 




R 


Non-IRDC 


G053. 11+00.05 


MM4 


19:29:20.4 


17:55:04 


21.7 


1.8 


40.0 


A 


IRDC 


G053. 11+00.05 


MM5 


19:29:26.3 


17:54:53 


21.7 


1.8 


45.0 


R 


IRDC 


G053. 25+00.04 


MM1 


19:29:39.0 


18:01:42 


24.9 


2.0 


35.0 


1 


IRDC 




IV J. 1V1 £l 


19:29:33.0 


18:01:00 


24.3 


2 
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G053. 25+00.04 


MM3 


19:29:44.0 


17:58:47 


23.5 


1.9 


35.0 


1 


IRDC 


G053. 25+00.04 


MM4 


19:29:34.5 


18:01:39 


24.5 


2.0 


37.0 


A 


IRDC 


G053. 25+00.04 


MM5 


19:29:39.4 


17:58:40 


23.7 


1.9 


35.0 


1 


IRDC 


G053. 25+00.04 


MM6 


19:29:31.5 


17:59:50 


23.8 


1.9 


46.0 


A 


IRDC 


G053. 31+00.00 


MM1 


19:29:50.0 


18:05:07 


22.0 


1.8 




i 


Blended 


G053. 31+00.00 


MM2 


19:29:42.1 


18:03:57 


25.2 


2.0 


25.0 


Q 


IRDC 


G053. 31+00. 00 


MM3 


19:29:49.7 


18:04:39 


22.2 


1.8 




i 


Blended 



of spectra for one clump of each proposed evolution- 
ary stage: the quiescent G028.53 MM3, the intermediate 
G025.04 MM4, the active G034.43 MM1, and the red 
G034.77 MM1, respectively. The Spitzer/IRAG image of 
the corresponding host IRDC is also displayed in each fig- 
ure. The active clump G034.43 MM1 (Figure g]) exhibits 
more molecular line emission than the quiescent clump 
G028.53 MM3 (Figure[2|). This pattern is typical through 
the sample. It is expected that when star-forming clumps 
evolve to later stages, they should display a rich molecu- 
lar spectrum with numerous lines because the proto-stars 
heat their environment releasing more complex molecules 
from dust mantles into the gas phase; in contrast, in ear- 
lier stages clumps should show a sparse spectrum with 
a smaller number of detectable molecular lines because 
many molecules are frozen out onto grain surfaces. How- 
ever, one must be aware that although the analysis of 
detection rates may reflect the chemical composition of 
clumps to some extent, it does not give unambiguous 
conclusions because detection rates may just show the 
dependence on H2 column density. A more definitive 
conclusion must be drawn from molecular abundances 
(see Sections 14.21 and |4.4[) . Figure [5] presents histograms 
which compare the number of detected molecular tran- 
sitions for each group. More molecular transitions are 
detected toward clumps that show signs of star forma- 
tion. The median number of detections for quiescent 
clumps is 6 molecular transitions. On the other hand, 
intermediate and active clumps have a median of 7 and 
9 detected molecular transitions, respectively. To quan- 
tify the differences between quiescent and active clumps, 
we have calculated the probability that these two dis- 
tributions in Figure |5] arise from the same underlying 
population using the Kolmogorov-Smirnov (K-S) test. 
The K-S test is a more robust method for measuring the 
similarity between two distributions than the comparison 
of the median and/or mean values. Through this paper, 
we will mention that the difference between two distribu- 
tions is statistically significant when the K-S test gives 
a percentage lower than 5%. The probability that the 



Note. — Units of right ascension arc hours, minutes and seconds, 
and units of declination arc degrees, arcminutes and arcscconds. 

a Velocity from the N 2 H+ JF t F = 123 012 transition. 

k Denotes the classification based on Spitzer/IRAC colors and the 

presence or absence of Spitzer/MlPS 24 )im emission. 

C Velocity and distance determined from the HCO^ line. 

^ Velocity and distance from fRathborne et al.l 1 120101) . 

c IRDC: the clump is located in an IRDC (the 92 sources of this study). 

Non-IRDC: the clump is not associated with an IRDC. Non-Detection: 

no molecular line was detected in this survey. Bright Blue: objects 

with bright 3.6 /im emission that arc predominantly uncxtinctcd stars. 

Blended: the angular separation between two sources is less than one 

Mopra beam. 

i i i i i 

Total Number of IRDC Clumps: 92 




HNC N 2 H + HC0 + HCN C 2 H H IJ C0 + HN 13 C HCjN HNCO SiO 
Molecular Lines 



Figure 1. Detection rates of the observed molecular lines toward 
IRDC clumps. The N2H + and C2H detection rates correspond to 
those of the JF 1 F = 123 - 012 and NJF = 1§2->0§1 transi- 
tions, respectively. Error bars are the Poisson noise (i.e., the root 
square of the number of detections). Values on bars are number of 
detections for a given molecule. 

Figures [2] [3l |H and [5] show representative examples 



Chemistry in Infrared Dark Cloud Clumps 7 
quiescent and active clumps are derived from the same parent population is 0.1%. 



Table 3 

N 2 H+ JFiF = 123 -> 012 Gaussian Parameters. 



IRDC 


Clump 


± rms 




Vl sr 


AV 


Name 


Number 


(K) 


(K) 


(km s — 1 ) 


(km s _1 ) 


G015. 05+00. 07 


MM1 


0.04 


0.38 ± 0.03 


25.21 ± 0.10 


3.18 ± 0.31 


G015.05+00.07 


MM2 


0.04 


0.19 ± 0.03 


29.87 ± 0.20 


3.64 ± 0.62 


G015.05+00.07 


MM3 


0.04 


0.18 ± 0.03 


28.41 ± 0.21 


2.50 ± 0.54 


G015.05+00.07 


MM4 


0.04 


0.44 ± 0.03 


27.95 ± 0.08 


2.37 ± 0.21 


G015.05+00.07 


MM5 


0.05 


0.25 ± 0.05 


29.93 ± 0.16 


2.03 ± 0.46 



Note. — This table is available in its entirety in a machine-readable 
form in the online journal. A portion is shown here for guidance re- 
garding its form and content. 

The median for red clumps is 6 detections. Red clumps 
(which show bright 8 /im emission) likely correspond to 
embedded Hn regions, a more evolved stage than active 
clumps. Despite this, they have the same median value 
as quiescent clumps. In the histogram, red clumps seem 
to have a bimodal distribution with a first population of 
8 clumps showing a low number of detected molecular 
transitions (left side of the histogram) and a second one 
of 7 clumps presenting a higher number (right side of the 
histogram). Moreover, we note that the first population 
only has one maser detection, and the second popula- 
tion has five ([Chambers ct al. 200§). These two different 
populations for red clumps may be the result of the evo- 
lutionary state of the H n region. We expect to detect 
more molecular lines from an early H n region compared 
to a late Hii region because, since the early Hn region 
has not had as much time to ionize the surrounding gas, 
it should have more molecular gas near the central heat- 
ing star(s). The molecular spectrum of a late stage Hn 
region, where much of the gas has been ionized, should 
be more sparse. 

The active clump G034.43 MM1 (Figure gj) exhibits 
a spectrum of a "hot core." It shows emission in the 
HNCO, HC 3 N, and CH 3 CN lines, which require high 
densities and temperatures for excitation (n cr it > 5 x 10 5 
cm~ 3 and E u /k > 10 K), and depend on the release of 
their parent molecules from dust grains to be formed. 
In addition, it shows SiO emission, which is usually as- 
sociated with molecular outflows si nce SiO abundance 
is highly enhanced by shocks (e.g., iSchilke et al.l 119971 : 
iCaselli et al.l 1199 7). In contrast, the quiescent clump 
G028.53 MM3 (Figure HJ) presents no emission of HNCO, 
HC 3 N, CH 3 CN, and SiO over the 3a limit. The inter- 
mediate clump G025.04 MM4 (Figure [3]) seems to be 
in an intermediate chemical state, only exhibiting one 
high-excitation line (HNCO). To carry out a more com- 
prehensive study of the full sample of IRDC clumps, we 
have split the histogram of detection rates (Figure [T]) into 
four new histograms. In Figure [7] the detection rates for 
each proposed evolutionary sequence are presented. The 
detection rates for the quiescent clumps are displayed in 
every panel for comparison. Quiescent, intermediate, ac- 
tive, and red clumps have 28, 23, 26, and 15 members, 
respectively. From the histograms in Figure it is seen 
that the HNC J = 1 ->• 0, N 2 H+ JF X F = 123 -> 012, 
and HCO + J = 1 — > lines are present in almost ev- 



ery IRDC clump at any evolutionary stage. Red and 
active clumps have a high detection rate (~95%) for the 
HCN J = 1 — > line, comparable with the three lines 
previously mentioned, while in less evolved clumps this 
value is lower (~69%). There is a high enhancement 
in the detection of the brightest C2H line, transition 
NJF — 1 1 2 — »• 1 1, toward active and intermediate 
clumps with respect to quiescent clumps. The percent- 
ages of detection for the H 13 CO + J = 1 —> line are 
roughly the same for red, intermediate, and quiescent 
clumps; however, the value for active clumps increases 
by ^30% with respect to the other three evolutionary 
stages. The detection rates of the HN 13 C J = 1 — > 
line show no significant difference between red and in- 
termediate clumps with respect to quiescent clumps; on 
the other hand, the value for active clumps increases by 
^30%. The most remarkable difference between detec- 
tion rates for the different evolutionary stages is found 
for the HC 3 N J = 10 -> 9 line. The detection of this 
line is significantly high for all stages with IR signs of 
star formation. HC3N is only detected in one quiescent 
source, and the difference between the detection rate in 
active and quiescent clumps is ^40%. While the HNCO 
Ji< a .K b — 4o,4 — > 3o,3 line is only detected twice (7%) in 
quiescent clumps, active and intermediate clumps have a 
HNCO JK a ,K b — 4o,4 — > 3o,3 detection rate of ~27%, in- 
dicating that this line is also observed more often toward 
clumps with current star formation. In the quiescent, in- 
termediate, and red clumps, SiO J — 2 — 1 emission was 
rarely detected. In active clumps, the detection rate of 
SiO reaches the 19%. 

3.3. Velocity Widths 

The widths of molecular lines give information about 
the turbulence of clumps. Although we cannot rule out 
that clumps have organized velocity structures, such as 
rotation, in this paper we will assume that turbulence 
dominates the widths of the lines. To investigate if tur- 
bulence is different in clumps with and without signs of 
star formation, we made histograms of the number distri- 
butions of full width at half maximum line widths (AV) 
for the four evolutionary stages using the three most com- 
monly detected lines: N 2 H+, HCO+, and HNC. In Fig- 
ured! the number distributions of N2H 4 " AV for each evo- 
lutionary sequence is presented. The range of line widths 
varies between 1.6 and 4.6 km s _1 for the full sample of 
IRDC clumps. The median values for quiescent, inter- 
mediate, active, and red clumps are 2.7, 2.8, 3.0, and 3.4 
km s _1 , respectively. Using the K-S test, we calculated 
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Figure 2. Quiescent clump example. Top: IRAC 3-color (3.6 fim in blue, 4.5 in green and 8.0 /an in red) image of the IRDC 
G028. 53-00. 25 overlaid with IRAM 1.2 mm continuum emissiom from Rathbornc ct al. (2006). Bottom: Molecular spectra of the quiescent 
MM3 clump. Quiescent clumps contain none IR emission. The circle shows the telescope beam size and marks the position of the observed 
clump. 
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Figure 3. Intermediate clump example. Top: IRAC 3-color (3.6 Atm in blue, 4.5 )im i n gre en and 8.0 fim in red) image of the 
IRDC G025. 04-00. 20 overlaid with IRAM 1.2 mm continuum emissiom from Rathborne ct al. (2006). Bottom: Molecular spectra of the 
intermediate MM4 clump. Intermediate clumps contain either a green fuzzy or a 24 /im source, but not both. The circle shows the telescope 
beam size and marks the position of the observed clump. 
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Figure 4. Active clump example. Top: IRAC 3-color (3.6 fim in blue, 4.5 lira in gree n and 8.0 )im in red) image of the IRDC 
G034. 43+00. 24 overlaid with IRAM 1.2 mm continuum emissiom from [Rathbornc ct al. (2006). Bottom: Molecular spectra of the active 
MM1 clump. Active clumps are associated with enhanced 4.5 (im emission, the so-called green fuzzies, and an embedded 24 ^m source. 
The circle shows the telescope beam size and marks the position of the observed clump. 
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Figure 5. Red clump example. Top: IRAC 3-color (3.6 am in blue, 4.5 )im in g reen and 8.0 fim in red) image of the IRDC G034.77- 
00.55 overlaid with IRAM 1.2 mm continuum emissiom from Rathbornc ct al. (2006). Bottom: Molecular spectra of the red MM1 clump. 
Red clumps are associated with bright 8 fim emission. The circle shows the telescope beam size and marks the position of the observed 
clump. 
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Quiescent Clumps: 28 
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Figure 6. Histograms with the number of detected molecular 
transitions per source for each evolutionary stage. The name of 
the corresponding evolutionary stage and the number of sources 
contained in each histogram are given on the top left side of each 
panel. We find, as expected, that more evolved clumps show a 
rich spectrum with numerous transition lines, except that the most 
evolved/red clumps show a population with few detections. 

a probability of 67% that active and quiescent, and 6% 
that red and quiescent distributions originate from the 
same parent populations. Although the K-S test shows 
that there is little difference between the active and qui- 
escent clumps, it seems that there is an increase in AV 
from clumps with no apparent star formation to active 
and red clumps. We think this trend is less clearly shown 
by the K-S test for active and quiescent distributions be- 
cause of the small number of members of each group; in 
spite of this, the K-S test shows a difference between the 
red and quiescent populations because the shape of the 
number distributions are more dissimilar. We conclude 
that the slightly broader widths found in the N2jH + line 
toward red and active clumps, compared with quiescent 
and intermediate clumps, is produced by the rise in the 
turbulence from the ongoing star formation activity. 

The number distributions of HCO + and HNC AV show 
no clear trend, and their histograms are not presented. 
The range of line widths varies between 2.0 and 13.0 
km s _1 for HCO + , and between 1.0 and 7.6 km s" 1 for 
HNC. The AV median values of HCO + for quiescent, 
intermediate, active, and red clumps are 3.8, 4.4, 3.6, 
and 4.0 km s _1 ; and for HNC are 3.2, 3.8, 3.4, and 2.8 
km s , respectively. We identify three reasons that may 
explain why HCO + and HNC do not present a similar 
trend to that one shown by N2H+. First, the HCO + and 
HNC lines are optically thick (see Section [4. 1 . 1[) and con- 
seque ntly, the line width s can be broadened by opacity 
(e.g.. iBeltran et al.l [20051) . Second, the HCO+ line can 
show broadening and wing emission produced by molec- 
ular outflow activity. Certainly, this is what is happening 
in G018.82 MM3, which has a AV of 13.0 km s" 1 in the 
HCO + line. Finally, because the Mopra beam is ^38" 
(0.8 pc at the average distance of 4.3 kpc), we are observ- 
ing the bulk motions of the gas and some level of clump- 
ing may also be present. This also could explain why the 



trend of broader N2H" 1 " lines toward more evolved objects 
is not as clear as we expect. However, si nce N2H" 1 " is 
thought to be a higher density tracer (e.g.. iCaselli et all 
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Figure 7. Detection rates of the observed molecular lines to- 
ward red, active, and intermediate clumps compared to detec- 
tions rates percentages in quiescent clumps. The N2jH + and 
C2H detection rates correspond to the JF\F = 123 — 012 and 
NJF = 1 2 2 -> 2 ltransitions, respectively. The evolutionary 
sequence name and the total number of sources are given on the 
top right corner of each panel. HNC and N2H + lines are detected 
in almost every IRDC clumps at every evolutionary stage. On the 
other hand, HC3N, HNCO, and SiO lines are predominantly de- 
tected in later stages of evolution. Error bars are the Poisson noise 
(i.e., the root square of the number of detections). 



I2002t iPirogov et al.l 12003ft , its emission should emanate 
mostly from the center of the clumps and thus we still 
can see a weak trend. 
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Table 4 

Detection Rates 







Quiescent 


Intermediate 


Active 


Red 


All 


Molecule 


Transition 


Clumps 


Clumps 


Clumps 


Clumps 


Clumps 



HNC 




J 











27 ± 5 (96%) 


22 ± 5 (96%) 


26 ± 5 (100%) 


15 ± 4 (100%) 


90 ± 10 (98%) 


N2H+ 


J 


= 1 — Q, F\ 


= 1 - 


l,F = 2 


— 2 




22 ± 5 (79%) 


22 ± 5 (96%) 


25 ± 5 (96%) 


9 ± 3 (60%) 


78 ± 9 (85%) 




J 


= 1 — 0, Fi 


= 2- 


1,F = 3 


— 2 




27 ± 5 (967o) 


26 ± 5 (1007o) 


2b ± 5 (1007o) 


13 ± 4 (877o) 


89 ± 9 (977o) 




J 


= 1 - 0, Fx 


= - 


l,F = 1 


— 2 




± vi (29%) 


12 ± 4 (52%) 


17 ± 4 (b57o) 


5 ± 2 (6670) 


42 ± 7 (467o) 


HCO+ 




J 


= 1 - 









or - I r* /riAfiy \ 

25 ± 5 (89%) 


17 ± 4 (74%) 


25± 5 (96%) 


14 ± 4 (93%) 


81 ± 9 (88%) 


HCN 




J 


= 1 - 









19 ± 4 (08%) 


Id ± 4 (707o) 


25± 5 (9d7o) 


14 ± 4 (93to) 


74 ± 9 (80to ) 


v. 2 [ ] 


N 


1 n 7 — 


3/2 - 


1/2, F = 


2 _ 


i 




j-O ^ L ±. yvo /o ) 




o ^ o ^uo /o y 


r.1 _|_ 7 ^crO^,^ 




V = 


1-0, J = 


3/2 - 


1/2, F = 


1 - 





1 ± 1 (4%) 


3 ± 2 (13%) 


10 ± 3 (39%) 


4 ± 2 (27%) 


18 ± 4 (20%) 




N = 


1-0, J = 


1/2- 


1/2, F = 


1 - 


1 


1 ± 1 (4%) 


6 ± 3 (26%) 


9 ± 3 (35%) 


3 ± 2 (20%) 


19 ± 4 (21%) 


H 13 CO+ 


N = 


1-0, J = 


1/2- 


1/2, F = 


0- 


1 


1 ± 1 (4%) 


1 ± 1 (4%) 


4 ± 2 (15%) 


2 ± 1 (13%) 


8 ± 3 (9%) 




J 


= 1 - 









13 ± 4 (46%) 


11 ± 3 (48%) 


20 ± 5 (77%) 


6 ± 3 (40%) 


50 ± 7 (54%) 


HN 13 C 




J 


= 1 - 









7 ± 3 (25%) 


8 ± 3 (35%) 


14 ± 4 (54%) 


5 ± 2 (33%) 


34 ± 6 (37%) 


HC 3 N 




J - 


= 10- 


- 9 






1 ± 1 (4%) 


4 ± 2 (17%) 


12 ± 4 (46%) 


5 ± 2 (33%) 


22 ± 5 (24%) 


HNCO 




J K a ,K b 


= 4(5,4 — 3o,3 






2 ± 1 (7%) 


6 ± 3 (26%) 


7 ± 3 (27%) 


3 ± 2 (20%) 


18 ± 4 (20%) 


SiO 




J 


= 2- 


1 






1 ± 1 (4%) 


± (0%) 


5 ± 2 (19%) 


2 ± 1 (13%) 


8 ± 3 (9%) 


H 13 CN 




J 


= 1 - 









± (0%) 


1 ± 1 (4%) 


3 ± 2 (12%) 


2 ± 1 (13%) 


6 ± 2 (7%) 


NH 2 D 




JK a ,K c 


= li,i — lo.i 






± (0%) 


± (0%) 


3 ± 2 (12%) 


± (0%) 


3 ± 2 (3%) 


SO 




Nj 


= 2 2 - 


- li 






± (0%) 


± (0%) 


± (0%) 


1 ± 1 (7%) 


1 ± 1 (1%) 


13 cs 




.1 


= 2- 


1 






± (0%) 


± (0%) 


2 ± 1 (8%) 


± (0%) 


2 ± 1 (2%) 


CH 3 CN 




J 


= 5 - 


4 






± (0%) 


± (0%) 


1 ± 1 (4%) 


± (0%) 


1 ± 1 (1%) 



_ Active CljrTpa 



_ Quiescent Clumps 



Figure 8. Histograms of the number distributions of N 2 H + 
linewidths (AV) for each evolutionary sequence. The name of the 
evolutionary stage is given on the top left corner of each panel. 
The vertical dashed lines indicate the median values of the AV for 
each distribution. Although the K-S test shows that active and 
quiescent clumps probably originate from indistinguishable parent 
populations, the K-S test shows that red and quiescent distribu- 
tions are different. The median values of the AV increase with the 
evolution of the clumps. 



Variations in N2H + line widths for different evolu- 
tionary stages towar d IRDC clumps ha ve also been ob- 
served previously by iSakai et al.l (120081) an d recently by 
iVasvunina et all (|2011[ ). ISakai et al.l (|2008l) observed 22 
MSX dark objects and 7 MSX sources located at dis- 
ta nces less tha n 4.5 k pc. MSX dark objects are defined 
by ISakai et al.l (|2008|) as objects with 8 /^m extinction 
features, although some of them have Spitzer 24 /xm 
pointlike sources; thus, MSX dark objects correspond 
to our intermediate and quiescent c lumps. The MSX 
sources correspond to our red clumps. ISakai et al.l (|2008l ) 
found that the AV of N2H" 1 " lines range between 0.6 and 



4.0 km s _1 , and the mean values for MSX dark objects 
and MSX sources are 2.3 and 2.5 km s , respectively. 
IVasvunina et al.l (|2011D observed 13 "quiescent" clumps 
(as defined by us), 11 "middle" clumps (which corre- 
spond to our intermediate and active clumps), and 13 
"active" clumps (which correspond to our red clumps). 
They found that the range of AV of N2jH + lines varies 
between 0.6 and 2.8 km s , and "quiescent," "middle," 
and "active" clumps have AV mean values of 1.4, 1.7, 
and 2.2 km s _1 , respectively. If the assumption that 
turbulence increases with star formation activity produc- 
ing broader molecular line widths in more evolved stages 
is correct, then the r e sults on N2H + line wid ths found 
bv ISakai et ail (|2008l ). IVasvunina et al.l (|201lD . and this 
work support the evo lutionary sequence proposed by 
IChambers eFail (pOOl ) for IRDC clumps. 

Figure |9] shows the AV of the most frequently de- 
tected molecular lines against the widths of N2jH + . The 
linewidths for HNC, HCO + , and SiO lines show a high 
degree of scatter. In the case of HNC and HCO + , this 
scatter is likely produced by the high line opacities (see 
Section [4.1 . 1|) and because these lines may also be trac- 
ing more external layers of the clumps. In the case of 
SiO, and also HCO + , the scatter is probably produced 
by molecular outflow activity, which increases the AV of 
the lines and adds extended wing emission. C2H appears 
to present the best correlation with N2H + . H 13 CO + , 
HN 13 C, and HC3N have slightly narrower line widths 
than N2H + , indicating that they trace similar optically 
thin gas emanating from the center of the clumps. 

From Figure |H1 it can be seen that there are four 
sources with very large linewidths in the HCO + (>8.0 
km s" 1 ) and SiO (>6.0 km s" 1 ) lines. They are G018.82 
MM3, G018.82 MM4, G035.59 MM1, and G033.69 MM2 
for HCO+ and G018.82 MM3, G019.27 MM2, G022.35 
MM1, and G028.04 MM1 for SiO. Notably, the four 
sources that have l arge SiO linewidths are active clumps. 
ISakai et al.l (|2010l ) mention that larger SiO linewidths 
and abundances are expected in younger sources than 
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Figure 9. Correlation plots of the N2H+ velocity widths against tJ 
lines are given on the top left corner of each panel. The dashed 1 
asterisks to intermediate clumps, green triangles to active clumps, a 
a high degree of scatter. C2H has the best correlation with N2H+ 
N 2 H+. 

MSX sources (our red clumps). However, as they also 
mention, strong shocks are necessary to produce SiO in 
the gas phase; therefore, some star formation activity is 
needed and the sources should be older than quiescent 
clumps. 

3.4. Kinematic Distances 

Because the N 2 H+ JFxF = 123 -> 012 line has a de- 
tection rate of 97% and shows no self-absorbed profiles, 
we used the N 2 H + systemic velocities of the 155 de- 
tected clumps from the initial sample to estimate kine- 
matic distances . The distances were obtained using the 
IClemensl (|1985f ) rotation curve of the Milky Way, scaled 
to (i?o, V ) = (8.5 kpc, 220 km s" 1 ) and included a small 
velocity correction (7 km s _1 ) accounting for the mea- 
sured solar peculiar motion. The distances of all clumps 
are listed in Table [5J To resolve the ambiguity be- 
tween the near and far kinematic distances, we used 
the fact that IRDCs are seen as dark extinction fea- 
tures against the bright mid-infrared background emis- 
sion. Hence, the ambiguity is resolved by assuming that 
all IRDCs are located at the near kinematic distance 



: velocity widths of several molecular lines. The name of the molecular 
le indicates unity. Blue squares correspond quiescent clumps, yellow 
d red circles to red clumps. Linewidths of HNC, HCO + and SiO show 
H 13 CO+, HN 13 C and HC3N have slightly narrower line widths than 

(e.g.. iSimon et~aTll2006bl Uackson et alj|200l . The dis- 
tances to the I RDC clumps in this work were previously 
established by Simon et al.l (|2006bl ). They carried out a 
morphological matching of the MSX mid-infrared extinc- 
tion with the 13 CO J = 1 — > molecular line emission 
from the Boston University-Five College Radio Astron- 
omy Observatory (BU-FCRAO) Galactic Ring Survey 
(GRS; Uackson et al.ll2006h . However, because the 13 CO 
J = 1 — > emission is a low density tracer (n cr i t = 2x 10 3 
cm ), there are typically several velocity components 
along one line of sight. Of the initial 155 detected IRDC 
clumps, the kinematic distance indicated by the high 
density tracer N 2 H+ (n crit = 3xl0 5 cm' 3 ) of 35 (23%) 
clumps differs from the distance inferred from 13 CO by 
ISimon et all (|2006b| ). Of these sources, 17 show strong 
signatures of star formation and are classified as red and 
active clumps. Because they are isolated clumps located 
apart from the IRDC and at a different distance, they 
are not physically associated with the IRDC and were 
excluded from further analysis in this paper. 
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4. DISCUSSION 



4.1. Derivation of Physical Parameters 

Assuming local thermodynamic equilibrium (LTE) 
conditions, the total column d ensity, N, can be derived 
from (e.g. JGarden et al.lll99i1) 



TV = 



8tw 3 Q rc 



exp(£ ( /fcT ex ) 



c 3 R g u A ui [1 - exp(-/ji//fcT ex )] 



t di; 



(1) 



where tv is the optical depth of the line, g u is the statis- 
tical weight of the upper level, A u \ is the Einstein coef- 
ficient for spontaneous emission, E l is the energy of the 
lower state, Q r ot is the partition function, v is the tran- 
sition frequency, and R is the relative intensity of the 
brightest hyperfine transition with respect to the others. 
R is only relevant for hyperfine transitions because it 
takes into account the satellite lines correcting by their 
relative opacities. It is 5/9 for N2H" 1 " and 5/12 for C2H. 
It is equal to 1.0 for transitions without hyperfine struc- 
ture. 

A particular case of equation [T] is the total column 
density of a linear, rigid ro tor molecule, which is given 
by (e.g.. iGarden et aJj fl991) 



N = 



3fc 



(T ox + hB/3k) exp(Ej /kT e , 



8n«B^R (J +1) [1 - exp(-hv/kT eK )] 



Tv dv , 

where /u is the permanent dipole moment of the 
molecule, J is the rotational quantum number of the 
lower state, Ej = hBJ(J + 1) is the energy in the level 
J, and B is the rotational constant of the molecule. 

To calculate column densities it is necessary to esti- 
mate the optical depths of the lines. We calculate the 
optical depths of N2H + and C2H using their hyperfine 
structure, and the optical depths of HCO + and HNC 
using their isotopologues, H 13 CO + and HN 13 C Under 
the assumption of LTE, all levels are populated accord- 
ing to the same excitation temperature (T ex ). We as- 
sume that T Py is equal to the dus t temperature (To) 
obtained bv iRathborne et a l. (2010). They determined 
Td for 59 sources and gave lower and upper limits for 
23 other sources in our sample of 92 IRDC clumps by 
fitting a graybody function to the spectral energy dis- 
tribution (SED). For the sources with lower and upper 
limits, the average value of the limits was used as Td- 
For the remaining 10 sour ces, without SEDs, we use 
the median values found bv IRathborne et all (|2010h for 
each evolutionary sequence: 23.7, 30.4, 34.5, and 40.4 K 
for quiescent, intermediate, active, and red clumps, re- 
spectively. When the weakest line needed to compute 
the optical depth was not detected above the 3<r level 
(N 2 H+ JF X F = 112 -> 012, C 2 H NJF = 1 § 1 -> \ 0, 
H 13 CO + , or HN 13 C), the intensity of the line was as- 
sumed to be 3T rms . This gives an upper limit of the line 
intensity and thus an upper limit in the optical depth. 

Actual optical depths were determined as follows. The 
main beam brightness temperature (T m b) is the antenna 
temperature (Ta) corrected by the main bean efficiency 
(Vmh), ?mb = TA/rjmb- The main beam brightness tem- 
perature of a line is related to the excitation temperature 
and the optical depth by 



T ulh = f[J(T ex )-J(T hg )](l-e 



(3) 



where / is the filling factor, t„ is the optical depth of 



the line, 7b g is the background temperature, and J(T) is 
defined by 

hv 1 



J(T) = 



k e hv / kT - 1 



(4) 



Taking the ratio of Equation [3] evaluated for two hy- 
perfine components or two isotopologues and assuming 
that the filling factors are the same, we obtain 



T, 



mbi 



1 



mb2 



J 2 (T ex ) - J 2 (T hg ) 



(5) 



where the subscripts "1" and "2" refer to two different 
hyperfine components or isotopologue molecules. Defin- 
ing "r" as the ratio between the optical depths, r = 
T 2/ T i, and utilizing the fact that the term between brack- 
ets is ~1, Equation [S] becomes 



1 



,-r 2 /r 



T, 



nib 1 



l-e 



mb 2 



(6) 



Solving this Equation numerically, we obtain the op- 
tical depth for one of the components, and by using the 
relationship between opacities we can obtain the optical 
depths of the second component. 

4.1.1. Optical Depths 

N 2 H+ 

The opacities of the N2H + hyperfine components are 
obtained from the ratio between the observed main beam 
brightness temperatures of the two brightest transitions: 
JFxF = 123 -> 012 and JF X F = 112 -> 012. Rewriting 
Equation [5J we obtain 



1 



g 5 123-012 



f mbn 



-Mnbia 



(7) 



where the subscripts "123 - 012" and "112 - 012" refer 
to the N 2 H+ hyperfine components JF 1 F = 123 -> 012 
and JF\F — 112 — > 012, respectively. We have also 
used the opacity ratio between these two hyperfine lines, 
\, which depends only on the transition mo- 



' 112-012 
ments. 

The N2H + opacities were calculated for 63 IRDC 
clumps and the values for the main component are listed 
in Table [5l The main component of N2H + has opacities 
less than 1.0 in roughly 60% of the sample. The other 
40% of the sample is moderately optically thick, with a 
maximum opacity of 4.6. The median values for quies- 
cent, intermediate, active, and red clumps are 0.8, 0.8, 
0.5, and 1.8, respectively. 

C 2 H 

The optical depths of the C2H hyperfine compo- 
nents are obtained using the intensity ratio between the 



NJF 



1S2 



Oil and NJF 



1 § 1 -> \ tran- 



sitions. The opacity ratio between these two lines is 
r= Ii2=<n_ = 2 (ITucker et al.lll97l . Thus, the opacity of 

T ll-00 

the brightest component is given by 



1 



1 - e T i2-oi T mbl2 _oi 



(8) 
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where the subscripts "12 — 01" and "11 — 00" refer to the 
C2H hyperfine components NJF — 1 1 2 — > 0^1 and 
NJF = l|l->0i0, respectively. 

The C2H opacities were calculated for 43 IRDC clumps 
and the values for the brightest transition are listed in 
Table [5] The emission of the brightest C2H hyperfine 
component is mostly moderately optically thick, with 
86% of the optical depths ranging from 1.0 to 6.0. The 
median values of the optical depths for quiescent, inter- 
mediate, active, and red clumps are 2.9, 2.4, 2.2, and 1.4, 
respectively. 

HCO+ and H 13 CO+ 



Rewriting Equation [51 we obtain 

l-e- T i2A T mb (H 13 CCT 



1 - e 



T mb (HCO+) 



(9) 



where "r" is the ratio between the optical depths, r = 
TI2/T13, and the subscripts "12" and "13" refer to the 
HCO + and H 13 CO + isotopologues, respectively. A way 
to determine "r" is to take the ratio of the expression for 
the column density (see Equation [2]) evaluated for both 
isotopologues 



HCO+ 
H 13 CO H 
[1 - exp(-/u> 12 /fcT ex )] 



(kT cx /hB 13 + 1/3) eME Jl3 /kT c: 
(kT c JhB 12 + 1/3) exp(S Ji2 /fcr e; 



(10) 



[1 - exp(-/u/ 13 /fcT cx )] 

where [HCO + /H 13 CO + ] is the isotopic abundance ra- 
tio. The isotopic abun dance ratio of [ 12 C / 13 C] ranges 
from -20 to -70 fe.g.. iSavage et all 120021) . Assuming 
a constant [HCO + /H 13 CO + ] abundance ratio of 50 for 
all sources, we compute the optical depths of 62 IRDC 
clumps. The values for HCO + are summarized in Ta- 



ble [5J All sources have optically thick HCO + emission, 
with opacities ranging between 7 and 75. On the other 
hand, H 13 CO + is optically thin in most of the sources 
(90%). The median values of HCO + optical depths for 
quiescent, intermediate, active, and red clumps are 30, 
20, 16, and 19, respectively; and the values for H 13 CO + 
are 0.6, 0.4, 0.3, and 0.4, respectively. 

HNC and HN 13 C 

Rewriting Equation |B] and following the same proce- 
dure used for HCO + , we obtain 



1 - e - T ^/ r 
1 — e~ r i2 



T mb (HN 13 C) 
T mb (HNC) 



(11) 



where r is the ratio between the optical depths, r 
^12/^13, and is given by 



(kT ex /hB 13 + 1/3) exp(£; ji3 /fcr cx ) 
(kT cx /hB 12 + 1/3) exp{E,/k,T cx ) 



HNC 
HN 13 C 



[1 - exp(-/u/ 12 /fc!r cx )] 
[1 - exp(-/ii/ 13 /fcT cx )] 



(12) 



where the subscripts "12" and "13" refer to the HNC and 
HN 13 C isotopologues, respectively, and [HNC/HN 13 C] is 
the isotopic abundance ratio. 

Assuming an [HNC/HN 13 C] abundance ratio of 50, we 
compute the optical depths of 85 IRDC clumps. The 
values for HNC are summarized in Table [5j All sources 
have optically thick HNC emission, whereas HN 13 C is 
optically thin in most of the sources (95%). The me- 
dian values of HNC optical depths for quiescent, inter- 
mediate, active, and red clumps are 24, 22, 18, and 16, 
respectively; and for HN 13 C are 0.5, 0.4, 0.3, and 0.3, 
respectively. 



Table 5 

Optical Depths and Column Densities. 



IRDC Clump 



Optical Depth 



r N 2 H+ r HCO+ TINC TC 2 H 



Column Density (cm ) 



N2H" 1 
xl0 1: 



HCO+ 
xlO 14 



HNC 
XlO 14 



C 2 H 
xlO 14 



HC3N 
xlO 13 



HNCO 
xlO 13 



SiO 
xlO 12 



G015.05 MM1 


1.00 


29.9 


26.9 


2.8 


1.56 


0.12) 


2.05(0.72) 


2.27 


1.09) 


2.75(3.02) 


1.67(0.69) 


G015.05 MM2 


0.04 




59.4 




0.71 


'0.15) 


0.58(0.21) 


2.66 


1.91) 






G015.05 MM3 


0.60 




28.6 




0.36 


0.10) 


1.04 


0.51) 






G015.05 MM4 






16.1 


1.5 


1.04 


0.12) 


0.58(0.20) 


1.56 


'0.65) 


3.06(3.29) 




G015.05 MM5 


0.36 


35.0 


52.2 




0.52 


0.16) 


1.01(0.56) 


2.40 


[1.55) 






G015.31 MM2 


2.89 




69.0 




1.18 


;i.25) 


2.19 


[1.79) 






G015.31 MM3 




50.2 


30.3 




0.31 


0.08) 


0.85(0.61) 


0.51 


0.51) 






G015.31 MM5 




38.2 


53.0 






1.10(0.63) 


1.36 


'0.95) 






G018.82 MM2 


2.16 


44.6 


43.5 




2.24 


'2.25) 


3.20(1.94) 


3.45 


2.10) 






G018.82 MM3 


2.27 


37.7 


27.9 




6.07 


'6.15) 


10.10(5.24) 


6.71 


3.14) 






G018.82 MM4 


0.12 


75.0 


45.3 




1.06 


'0.18) 


8.50(6.72) 


5.02 


2.85) 




0.98(1.04) 


G018.82 MM6 


0.12 








1.22 


'0.12) 


0.52(0.19) 


0.80 


0.35) 






G019.27 MM2 


0.81 




33.8 


5.6 


2.18 


'0.15) 




4.09 


1.40) 


4.15(4.59) 




G022.35 MM1 


0.31 


14.4 


20.9 




0.40 


,0.10) 


0.85(0.33) 


1.13 


0.47) 






G022.35 MM2 






33.2 










2.02 


'1.07) 






G023.60 MM7 


1.64 




27.6 




4.13 


'4.07) 




3.50 


1.67) 






G023.60 MM9 


0.56 




29.3 


2 1 


0.68 


,0.11) 




1.47 


,0.71) 


1.51(2.16) 




G024.08 MM2 


0.99 


43.0 


26.4 




0.70 


'0.12) 


1.77(0.71) 


2.51 


'1.12) 




G024.08 MM3 


1.27 


31.5 


24.1 




0.93 


'0.92) 


1.78(0.89) 


1.91 


'0.85) 






G024.08 MM4 


0.27 


66.5 


31.1 




0.48 


0.13) 


2.43(2.20) 


1.84 


0.94) 






G024.33 MM2 






25.5 




3.06 


0.18) 


1.37(0.36) 


6.16 


1.17) 




4.27(1.10) 


G024.33 MM3 








1.9 


2.65 


;o.i7) 


1.38(0.35) 


0.85 


'0.29) 


3.41(4.60) 


2.81(0.90) 3.96(1.41) 


G024.33 MM4 










1.57 


'0.13) 


0.79(0.30) 


0.78 


'0.28) 






G024.33 MM5 


0.26 






3.2 


2.17 


0.15) 


2.81(0.46) 


0.79 


0.29) 


2.63(3.99) 


5.44(0.94) 
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IRDC Clump Optical Depth Column Density (cm 2 ) 

tn 2 h+ t hco+ thnc -t C2 h N 2 H+ HCO^ HNC C2H HC3N HNCO SiO~ 

xlO 13 xlO 14 xlO 14 xlO 14 xlO 13 xlO 13 xlO 12 



G024.33 


MM7 





67 


15.0 


10.3 


1.2 


G024.33 


MM8 






13.2 


15.8 




G024.33 


MM9 






10.9 


11.5 




G024.33 


MM11 





59 


29.9 


23.9 




G024.60 


MM2 





92 


22.4 


16.4 


1.1 


G025.04 


MM2 






13.2 


18.3 




G025.04 


MM4 





06 


8.7 


24.6 


1.0 


G027.75 


MM2 


3 


17 








G027.94 


MM1 


1 


39 




22.9 


3.1 


G028.04 


MM1 


1 


15 


24.1 


60.7 


1.3 


G028.08 


MM1 


1 


25 




40.7 




G028.23 


MM1 


3 


62 








G028.28 


MM4 





88 


13.9 


11.1 


3.0 


G028.37 


MM1 








20.6 


1.9 


G028.37 


MM2 





46 




24.6 


3.0 


G028.37 


MM4 





06 


12.8 


11.8 


4.8 


G028.37 


MM6 





19 


22.2 


20.3 


2.3 


G028.37 


MM9 






21.9 


15.0 




G028.37 


MM11 





23 




41.8 


3.3 


G028.37 


MM12 


1 


37 


65.8 


74.0 




G028.37 


MM13 


3 


32 








G028.53 


TV T TV T O 

MM3 


1 


73 


19.4 


23.7 




G028.53 


MM5 


4 


45 


46.0 


46.7 




G028.53 


MM7 








49.0 




G028.53 


TV T TV T O 

MM8 






69.4 


37.1 


1.5 


G028.53 


TV T TV /Tn 

MM9 





81 


13.9 


15.0 




/"i rvo r" 

G028.53 


TV J TV K 1 rv 

MM 10 





10 


13.7 


17.2 


2.9 


G028.67 


MM1 






31.0 


31.2 




G028.67 


MM2 





83 




47.7 




G030.14 


MM1 


1 


58 


37.8 


32.6 




/ A j \ • \ j \ f f-r 

G030.57 


MM1 





11 




19.6 


1.1 


G030.57 


TV J TV T O 

MM3 





42 




45.2 


6.0 


G030.97 


MM1 






16.4 


9.4 




G031.97 


MM5 





04 




15.2 


3.5 


G031.97 


MM7 





45 


7.3 


14.5 


3.5 


G031.97 


TV J TV T O 

MM8 








9.5 




G033.69 


MM1 






9.8 


11.4 




G033.69 


TV T TV TO 

MM2 






17.1 


14.6 


0.8 


G033.69 


MM3 






18.5 


16.4 




G033.69 


MM4 






12.4 


16.2 




G033.69 


MM5 
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4.1.2. Column Densities 

The column densities of N 2 H+, HCO+ and HNC are 
calculated using equation [2J The SiO molecule is also 
a linear, rigid rotor. Thus, we can also use equation [2] 
to estimate SiO column densities. However, since we 
observed only one SiO isotopologue, we need to assume 
that the SiO emission is optically thin and the filling 
factor is 1.0 (which gives us the beam-averaged column 
density). Therefore, replacing equation [3] in equation^ 
we obtain 



N = 



3fc (T cx + hB/3k) exp(Ej/kT ex ) 



8ir 3 Bn 2 R 
1 



[J(T ex ) - J(T bg )] 



(J + l) [1- 
T mh dv 



exp(-hv/kT cx )] 



(13) 



The values of the parameters used in the column den- 
sity calculations (permanent dipole moment, rotational 
constant) of N 2 H+, HCO+, HNC, and SiO are summa- 
rized in Table [6] 



Table 6 

Parameters used for N 2 H+, HCO+, HNC, and SiO Column Density 
Calculations 

Molecule Dipole Moment Rotational Constant References 
a (D) B (GHz) 



N 2 H+ 

HCO+ 

H 13 CO+ 

HNC 

HN 13 C 

SiO 



3.40 
3.89 
3.89 
3.05 
3.05 
3.10 



46.586871 
44.594423 
43.377302 
45.331980 
43.545600 
21.711979 



1, 2, 3 
4, 5, 6 
4, 5, 6 
7 ,8 
7 ,8 
9, 10 



Not e. 
fl990T) : 



(1) jBotschwin3 i 19841): (21 IHavenith et al.l 

(3) Pagani ct al. (2009); (4) Botschwina ot al. (1993); 



References. 



(5) | Yama guchi et af| 
Blackroa^^^^D ]197€ 



Raymonda ct al. (1970 



(1994) 

(8) 



(6) ILattanzi et al.l [2007); 



| van der Tak et al l (1200 
(10) IMoUaaghababa et al.l J199lT) 



(7) 
(9) 



The column densities of the remaining molecules (C2H, 
HC3N, HNCO) are determined using equation [TJ Al- 
though C2H is a linear molecule, its rotational energy 
levels are described by the rotational quantum number 
N, instead of J. HC3N is also linear and its partition 
function is the sam e used for diatomic, linear molecules 
(Blak e et al.l [l987h ; we have not included the factor 3 
due to the nuclear spin degeneracy because the hyperfine 
structure is not resolved. On the other hand, HNCO is a 
standard asym metric top molecu le with a different par- 
tition function ([Blake et al.lll987l ). For C 2 H the column 
density was calculated using the optical depths obtained 
in Section ElTLTl For HC 3 N and HNCO, the emission 
was assumed to be optically thin and the filling factor 
was used to be equal to 1.0. The values of the parame- 
ters used in the column density calculations (partition 
function, rotational constants, statistical weight, Ein- 
stein coefficient for spontaneous emission, and energy of 
the lower state) of C 2 H, HC3N and HNCO are summa- 
rized in Table [7] 

Upper limits to the column densities of the most com- 
monly detected molecular lines (N 2 H + , HCO + , HNC) 
were estimated when they were not detected over the 3cr 
level. In the case of the optically thick HCO + and HNC 
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Figure 10. Histograms of the number distributions of N2H+ col- 
umn densities for each evolutionary sequence. The name of the 
evolutionary stage is given on the top left corner of each panel. 
The vertical dashed lines indicate the median values of the column 
densities for each distribution. The diagonal patterns indicate the 
portions of the histograms that correspond to upper and lower lim- 
its (see bottom panel, right side). The median values and the K-S 
test suggest that column densities increase with the evolution of 
the clumps. Median values of the column densities are given in 
Tabled 

emission lines, the limits were obtained first for their iso- 
topologues (H 13 CO + and HN 13 C). Then by using the iso- 
topic abundance ratio we calculated the column densities 
for HCO + and HNC. In all non-detections, the emission 
was assumed to be optically thin, with T m b = 3T rms , and 
AV equal to 2.0 km s" 1 . 

Additionally, lower limits to the column densities were 
estimated for N 2 H + and HCO + when our normal ap- 
proach was not applicable. For N 2 H + , the intensity ra- 
tio between the two brightest components, which is used 
to derive opacities, should not be higher than its value 
in the optically thin limit, 5/3. However, in some cases, 
this ratio exceeds the optically thin limit, but neverthe- 
less remains consistent with 5/3 within the uncertain- 
ties. To calculate a lower limit for the column density, 
the main N 2 H + component was assumed to be optically 
thin. For HCO + , in several cases the spectra show self- 
absorbed profiles, and no Gaussian fit was performed. 
Lower limits to the column densities were computed us- 
ing the H 13 CO + emission line and assuming that its emis- 
sion is optically thin. Then, using the [HCO+/H 13 CO+] 
isotopic abundance ratio, we obtain the HCO + column 
densities. 

Derived column densities using the optical depths ob- 
tained in Section [4.1.11 were multiplied by the filling fac- 
tor to get the beam-averaged column density. The filling 
factor was calculated from Equation [3l Column densities 
for N 2 H+, HCO+, HNC, SiO, C 2 H, HC 3 N, and HNCO 
molecules (including all limits) are listed in the Table [5j 
Median values of the column densities for each evolu- 
tionary sequence (excluding limits when the lines are not 
detected above the 3a level) are shown in Table [5] 

N 2 H+ 

Figure [TU] shows the number distributions of the N 2 H + 
column densities for each evolutionary sequence. The K- 
S test gives a probability of 0.09% that quiescent and ac- 
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Figure 11. Histograms of the number distributions of HCO + 
column densities for each evolutionary sequence. The name of the 
evolutionary stage is given on the top left corner of each panel. 
The vertical dashed lines indicate the median values of the column 
densities for each distribution. The diagonal patterns indicate the 
portions of the histograms that correspond to upper and lower 
limits (see bottom panel, right side). The median values and the 
K-S test suggest that column densities increase with the evolution 
of the clumps. Median values of the column densities are given in 
Table [8] 

tive clump distributions originate from the same parent 
population. Figure [TU] shows a trend of increasing col- 
umn densities with the evolution of the clump, which is 
confirmed by the median values (see Table [3]) and the K- 
S test probability. This trend does not necessary imply 
a change in the chemistry. An increase of the N2H + col- 
umn density could simply indicate that the total H2 col- 
umn density increases while maintaining constant chem- 
ical composition. In Sections 14.21 and 14.41 we discuss 
the molecular abundances, which more definitively show 
chemical variations. This argument applies to all the 
molecules discussed in this section. Sa kai et ail (2008) 
obtained N2H + column densities over a sample of 11 
Spitzer and MSXdark objects (no 24 and 8 /im emission, 
respectiv ely) that we c a n com pare with our quiescent 
clumps. iPirogov et al.l (|2003f ) obtained N 2 H + column 
densities of 34 clumps, where massive stars have already 
formed, that we can compare with our active clumps. 
Table IH1 summarizes the results of previous works, which 
are in good agreement with our results. 

HCO+ and H 13 CO+ 

Figure UTI shows the number distributions of the HCO + 
column densities for each evolutionary stage. The K-S 
test gives us a probability of 2% that quiescent and ac- 
tive distributions originate from the same parent pop- 
ulation. Figure [TT] shows a trend of increasing column 
densities with the evolution of the clump, which is sup- 
ported by the me dian values ( s ee Ta ble [5t and the K- 
S test probability. ISakai et afl (|2010f ) and iPurcell et al.l 
(|2006h obtained H 13 CO+ and HCO+ column densities to- 
ward 20 and 79 massive clumps, respectively, associated 
with IR emission, methanol masers and UC H 11 regions. 
We compare this sample with our active and red cores 
in Table [9] The hig h HCO + column densities obtained 
by ISakai et al.l pOlOh and IPurcell et all ((2001I confirm 



Figure 12. Histograms of the number distributions of HNC col- 
umn densities for each evolutionary sequence. The name of the 
evolutionary stage is given on the top left corner of each panel. 
The vertical dashed lines indicate the median values of the column 
densities for each distribution. The diagonal patterns indicate the 
portions of the histograms that correspond to upper and lower 
limits (see bottom panel, right side). The median values and the 
K-S test suggest that column densities increase with the evolution 
of the clumps, from quiescent, to intermediate, to active clumps. 
Median values of the column densities are given in Table [8] 

the trend that HCO + column densities increase with the 
evolution of the clumps, as seen in Figure [TT] 

HNC and HN 13 C 

Figure [12] shows the number distributions of the HNC 
column densities for each evolutionary stage. The K-S 
test gives a probability of 5% that quiescent and active 
distributions originate from the same parent population. 
Figure [12] shows a weaker trend (compared with N2H+ 
and HCO + ) of increasing column densities from quies- 
cent to active clumps. There are only a few systematic 
studies of HNC toward a l arge sample of ma ssive star 
forming or IRDC clumps. ISakai et ail (|2010l ) observed 
HN 13 C toward 20 massive clumps that we can compare 
with our active clumps. Their values, showed in Table [HI 
are in good agreement with our results. 

C 2 H 

Figure [TJ|] shows the number distributions of the C2H 
column densities for each evolutionary sequence. The 
K-S test gives us a probability of 38% that active and 
quiescent distributions originate from the same parent 
population. Figure [T3] shows that the trend of increas- 
ing column densities with the evolution of the clumps 
is less evident for C2H than for previous molecules. Al- 
though the median values obtained for each evolutionary 
sequence support this trend, the K-S test probability 
is not low enough to be confident. Table [2] compares 
the column densit ies obtained for 20 massive clumps by 
ISakai et al.l (pOlOh with our active clumps, which are in 
a reasonable agreement. 

HC 3 N 

The median values for the HC3N column densities are 
summarized in Table [8] HC3N emission was detected 
only in one quiescent clump, G028.37 MM12, whose col- 
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Figure 13. Histograms of the number distributions of C2H col- 
umn densities for each evolutionary sequence. The name of the 
evolutionary stage is given on the top left corner of each panel. 
The vertical dashed lines indicate the median values of the column 
densities for each distribution. The diagonal patterns indicate the 
portions of the histograms that correspond to upper limits. The 
median values slightly increase with the evolution of the clumps; 
however, the K-S test shows this variation is not statistically signif- 
icant. Median values of the column densitie s are given in Table [8] 

umn density was 3.05 x 10 12 cm' 2 . ISakai et al.l (2008) 
obtained HC3N column densities for 7 Spitzer and MSX 
dark objects (no 24 and 8 /an emission, respectively). 
Their values are comparable to that found in our one 
quiescent clump (see Table [9]). 

HNCO 

The median values for t he HNCO column densi ties are 
summarized in Table HI iZinchenko et al.l (|2000h calcu- 
lated HNCO column densities toward 20 massive star 
forming clumps using rotational diagrams. These sources 
were selected based on the presence of water masers and 
strong SiO emission, indicating they are currently form- 
ing stars. They obtained higher values than those found 
in our IRDC sample (Tabled]). This may be explained 
by the war mer rotational temperatu res (median of 88 K) 
obtained by IZinchenko et al.l (|2000D . 

SiO 

The median values for the SiO column densities are 
summarized in Table [8] SiO emission was detected in no 
intermediate clumps and only in one quiescent clump, 
G024.33 MM7, wh o se co lumn density is 1.36 x 10 12 
cm- 2 . ISakai et al.l (|2010f ) and iMiettinen etaLl (|2006f ) 
obtained SiO column densities for 17 and 15 high-mass 
star-forming clumps, respectively, associated with IR and 
maser emission and UC H11 regions (See Table |9|). Their 
values agree with those we calculated for active clumps. 
Red clumps, however, have lower column densities. 

4.2. Chemistry 

In order to study the chemistry toward the IRDC 
clumps, and investig ate if the different evol utionary se- 
quences proposed by IChambers et all (|2009T ) are chemi- 
cally distinguishable, we calculate molecular abundances 
using the 1.2 mm co ntinuum emission obtained by 
iRathborne et al.l (|2006l ) and abundance ratios of selected 
species. 



4.2.1. Molecular Abundances 

To estimate molecular abundances with respect to 
molecular hydrogen, we take the ratio between the col- 
umn density of a given molecule and the H2 column 
density derived from dust emission. Observations of 
1.2 mm dust c o ntinuu m emission were obtained by 
IRathborne et al.l ((2006) with an angular resolution of 
11". In order to calculate the molecular abundances, 
we smoothed the continuum emission data to the angu- 
lar resolution of the molecular line data at 3 mm, 38". 
The resulting 1.2 mm peak flux, H2 column densities 
and molecular abundances for the N 2 H + , HCO + , HNC, 
C 2 H, HC 3 N, HNCO, and SiO lines (including limits for 
N 2 H+, HCO+, and HNC) are listed in Table [TO] In gen- 
eral, the number distributions of abundances show few 
differences between the evolutionary sequences, with the 
exception of N 2 H+ and HCO+ (Figure ITU and IT5) . The 
K-S test gives a probability of 0.09% and 1.2% that qui- 
escent and active distributions originate from the same 
parent populations for N 2 H + and HCO + , respectively. 
On the other hand, the evolutionary status of the clumps 
is not evidently distinguishable from the abundances of 
HNC, C 2 H, HC 3 N, HNCO, and SiO, although their me- 
dian values change with the evolution of the clumps. Me- 
dian values of the molecular abundances for each evolu- 
tionary sequence (excluding the limits when the lines are 
n ot detected) are shown in Table [SJ 

IVasvunina et al.l (|2011f) observed 37 IRDC clumps in 
the fourth quadrant with the Mopra telescope. They 
also obtained molecular abundances for N 2 H + , HCO + , 
HNC, SiO, C 2 H, HC 3 N, and HNCO. Despite a couple 
of different assumptions (jVasvunina et al.l (|2011[) did not 
smooth their 1.2 mm dust continuum emission to the an- 
gular resolution of Mopra (38") and they used NH 3 tem- 
pearatures), the abundances they find are in good agree- 
ment with the abundances determined in this work, ex- 
cept for HNC. They determined HNC abundances about 
an order of magnitud lower than our values. Because 
IVasvunina et al.1 (|2011h did not have the isotopologue 
HN 13 C, they assumed that the emission was optically 
thin in order to get an estimate of the column density. 
However, as we discussed in Section 14.1.11 the optically 
thin assumption is not valid for HNC in IRDCs because 
the emission is optically thick. 



4.2.2. Abundance Ratios 

Molecular abundance ratios that can be used to es- 
timate the age and mark the evolutionary stages of 
star-forming regions are known as "chemical clocks." 
Only molecules that show differential abundances with 
time can be used to evaluate the evolutionary sta- 
tus of a star-forming region. Chemical clocks have 
been studied in depth in low-m ass star-forming regions 
(e.g. jEmprechtinger et a l. 2009); the concept of chemical 
clocks has been less developed in the context of high-mass 
star-forming regions. Recent studies in small samples of 
high-mass star-forming regions show that chemical clocks 
can be extended t o their high-mass counterpart (e.g., 
iFontani et al.ll20Tll) . In this section, we explore the idea 
that N 2 H+/HCO+ and N 2 H+/HNC ratios can be used 
as chemical clocks. We choose this set of lines because 
they reveal the clearest significant chemical variations. 
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Table 7 

Parameters used for C2H, HC3N and HNCO Column Density Calculations 



Molecule 


Qrot 




Rot. Const. 
(GHz) 


9u 


A u i 
(xlO- 6 s" 1 ) 


EJk 


C 2 H 


kT ex /hB + 1/3 


B 


= 43.674518W 


5.0 


1.52757 


0.00216 


HC 3 N 


fcT ex //iB + 1/3 


B 


= 4.5490586< 2 ) 


21.0 


58.1300 


19.6484 


HNCO 


MfcTcx^/^ABC)] 1 / 2 


A 


= 918.417805( 3 ' 


9.0 


8.78011 


6.32957 




B 


= 11.071010< 3 ) 












C 


= 10. 910577( 3 ) 









Note. — Refe rence s. (1) IPadovani et aU <2009l) ; (2) 
IThorwirth et all 12000) ; f 31 ILapinov et al.l 1 120071 ). Values for g„, A u , 
and E L /k were obtain ed from The Cologne Da tabase for Molecular 
Spectroscopy (CDMS) IMiiller et alJI200U[2005l'l . 

N 2 H+ and HCO+ 

It has been suggested that the line and abundance 
ratios of N2H + and HCO + are anticorrelated, varying 
spatially over a given cloud and for different sources. 
Th e change in the N2H + /HC O + ratio was firs t note d 
bv lTurner fc Thaddeusl (fl977h an dlS nvder et al.1 (fl977[) . 
and recently by IKim et al.l (|2006f ) and ILo et all ((2009); 
however, it has not been studied in a large sample of 
star-forming clumps. Although we cannot examine spa- 
tial variations, due to our single-pointing observations, 
we can test whether the abundance ratios show dif- 
ferences for the four evolutionary stages proposed by 
iChambers et~aT1 (pOul ). 

According to ion-molecule schemes of dense molecula r 
clouds (e.g.. iSnvder et al.l H977I: Uorgensen et al.l l2004n . 
N2H + and HCO + are formed primary through reactions 
with H^": 



H+ 



CO -> HCO+ + H 2 



(14) 
(15) 



The main destruction mechanisms for N2H + are elec- 
tron recombination, when CO is frozen out on dust 
grains, 

N 2 H+ + e -> N 2 + H , (16) 

and by reacting with CO, when CO is found in the gas- 
phase at standard abundances QCO/H2] ~ 10~ 4 ), 



N 2 H+ + CO ->• HCO+ + N 2 



(17) 



The dominant removal mechanism for HCO + is elec- 
tron recombination 



HCO H 



CO + H . 



(18) 



Ch emical models of low-mass star-for ming regions 
(e.g. JLee et al.ll2004t Uorgensen et al.ll2004f) a nd one of a 
massive star-forming region ( Kim et al.ll2006f) show that 
the N 2 H + /HCO + ratio is high in early stages due to 
HCO + depletion (because most of the CO is frozen out) 
and is low in later stages when CO is evaporated from 
the dust grains (T ~20-25 K). Because CO is both the 
main supplier of HCO + and the main destroyer of N2H + , 
the theory suggests that the HCO + abundance increases 
with respect to N2H + as the clump evolves to a warmer 
phase. It is notable, for reasons not yet fully understood, 



that N2H + is resistant to depletion in cold, dense re- 
gions. For a long time, it was thought that N2, the main 
supplier of N2H+, had a lower binding strength to the 
grain surface than CO. However, recent measurements 
show that they are similar (jOberg et al.l I2005T ) ■ Hence, 
the chemical network that dominates the formation and 
destruction of N2H+ and HCO + in early stages of star 
formation is given by Equations Q31 HS1 HH and [T51 but 
for more evolved stages, Equation [TC] becomes less im- 
portant and is replaced by Equation 1171 A complete 
chemical network of Hjj~, N2H" 1 ", and HCO + , including 
for mation and destructio n rates, is presented in detail 
bv I Jorgensen et al.l (|2004D . 

Figure 1161 shows the number distributions of the 
N2H + /HCO + abundance (or column density) ratio for 
each evolutionary sequence. The median values for each 
distribution are 0.08, 0.12, 0.10, and 0.07 for quiescent, 
intermediate, active, and red clumps, respectively. As 
can be seen from the Figure 1161 quiescent clumps do not 
present the predicted abundance ratio to follow the trend 
of decreasing N2H + /HCO + ratio with the evolution of 
the clumps; their ratios are more randomized. For this 
reason, we use the K-S test for the intermediate and red 
clump distributions. The test shows that the probability 
of both populations being the same is low (3%). The val- 
ues of the median and the K-S test probability support 
the theoretical idea that the N2H + /HCO + abundance 
ratio acts as a chemical clock (at least for the interme- 
diate, active, and red clumps); however, why this trend 
does not extend to quiescent clumps is unclear. We note 
that the relative abundance between N2H + and HCO + 
acts as a chemical clock despite the abundances of both 
molecules independently increasing with the evolution of 
the clumps. We stress that a large number of values 
used in the histograms are limits and that a number of 
assumptions have gone into calculating the column den- 
sities and abundances. A sample with a larger number 
of strong detections would be very helpful for confirming 
these trends. 

The trend of decreasing N2H + /HCO + ratio with the 
evolution of the clumps may not extend to quiescent 
clumps due to one or a combination of the following rea- 
sons: (a) the sample of quiescent clumps could contain 
clumps with embedded star formation that was not de- 
tected by Spitzer, (b) CO has not had time to freeze 
out on dust grains, (c) the theoretical predictions at the 
very early stages are incorrect, and/or (d) due to the 
large beam of Mopra, we are probably detecting con- 
taminating, diffuse emission which is not directly related 
to the star-forming process in the center of the clumps. 
Considering the last point, since the size of the Mo- 
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Figure 14. Histograms of the number distributions of N2H+ 
abundances for each evolutionary sequence. The name of the evo- 
lutionary stage is given on the top left corner of each panel. The 
vertical dashed lines indicate the median values of the abundances 
for each distribution. The diagonal patterns indicate the portions 
of the histograms that correspond to upper and lower limits (see 
bottom panel, right side). The median values and the K-S test sug- 
gest that abunda nces increase w i th the evolution of the clumps, as 
it is also found by Busquct ct al. (2011) in their chemical modeling. 
Median values of the column densities are given in Table [8] 

pra beam is 38", we are typically observing emission 
on physical scales ~0.8 pc. Typical sizes of clumps in 
hig h-mass star-forming regions and in IRDCs are ^0.4 
pc (iFaundez et al.ll2004D and ~0.5 pc (|Rathborne et al.l 
2006), respectively. This shows that we are not only trac- 
ing the densest regions associated with star formation, 
but a lso the environment between clumps. ILee et aLl 
( 2004) predict that the variation of the molecular abun- 
dances is much larger when we consider small scales in- 
side the clumps. Further studies at higher angular res- 
olution, using interferometers, are needed to understand 
the discrepancy shown by quiescent clumps. 

N 2 H+ and HNC 

Figure [T7l shows the N2H+/HNC abundance ratio for 
each evolutionary stage. The median values for each dis- 
tribution are 0.06, 0.06, 0.08, and 0.08 for quiescent, in- 
termediate, active, and red clumps, respectively. The 
K-S test for the quiescent and active clump distributions 
shows that the probability of both populations being the 
same is low (5%). The median values and the K-S test 
suggest that there is a trend of increasing N2H + /HNC 
abundance ratio with the evolution of the clumps. We 
note that N2H+ and HCO + are strongly linked because 
of their pathways of formation and destruction; on the 
other hand, a connection between N2H 4 " and HNC is less 
evident. The trend between these two molecules suggests 
that HNC may be preferentially formed in cold gas. 

Other Abundance Ratios 

It is expected that abundances of SiO and more com- 
plex molecules, such as HC3N and HNCO, increase in 
more evolved star-forming regions. In evolved regions, 
these molecules or their parent molecules are released 
from the dust grains by shocks and sublimation due to 
the temperature increase. In fact, the abundance ratio 



Figure 15. Histograms of the number distributions of HCO + 
abundances for each evolutionary sequence. The name of the evo- 
lutionary stage is given on the top left corner of each panel. The 
vertical dashed lines indicate the median values of the abundances 
for each distribution. The diagonal patterns indicate the portions 
of the histograms that correspond to upper and lower limits (see 
bottom panel, right side). The median values and the K-S test 
suggest that abundances increase with the evolution of the clumps. 
Median values of the column densities are given in Table [8] 

between each one of them and N2H" 1 ", HCO + , or HNC 
show higher ratios in later stages of evolution. Because 
the number of SiO, HC3N and HNCO detections is low, 
the histograms with the number distributions of abun- 
dance ratios were not made. On the other hand, C2H 
has a larger number of detections and shows no clear 
trends when abundance ratios are taken with respect to 
N 2 H+, HCO+, or HNC. 

4.3. Assumptions and limitations in calculating column 
densities and abundances 

This section discusses the major assumptions and lim- 
itations that affect the calculations of column densities 
and molecular abundances. Although the assumptions 
used in this work are widely used in the literature, it is 
important to be aware of their implications. 

Temperature. Under the assumption of local ther- 
modynamic equilibrium (LTE) conditions, we assume 
that the excitation temperature for all molecules is iden- 
tical to the dust temperature for a given source. This as- 
sumption is reasonable since the critical densities of the 
molecules that were used to calculate physical parame- 
ters (between 2 x 10 5 - 2 x 10 6 cm" 3 ) are comparable to 
the d ensity of the sources (~10 6 cm~ 3 ; iRathborne et al.l 
I2010T ). implying they are close to thermalization. How- 
ever, the dust temperature is an average temperature 
along the line of sight, and may not necessarily reflect 
the temperature of the dust and gas in the portion of 
the clumps/cores which are dense enough to produce the 
emission lines we are observing. We note that the molec- 
ular abundance ratios (N 2 H+/HCO + and N 2 H+/HNC) 
are practically independent of temperature because of 
a similar dependence on temperature for the individual 
abundances. The abudance ratio for a given source at 15 
and 50 K varies only ~0.4%. 

Isotopic abundance ratio. In the calculation of op- 
tical depths and column densities for HCO + and HNC, 
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Figure 16. Histograms of the number distributions of 

N2H + /HCO + abundance ratios for each evolutionary sequence. 
The name of the evolutionary stage is given on the top left cor- 
ner of each panel. The vertical dashed lines indicate the median 
values of the abundance ratios for each distribution. The diagonal 
patterns indicate the portions of the histograms that correspond 
to upper and lower limits (see bottom panel, right side). The me- 
dian values and the K-S test suggest that N2H + /HCO + ratios act 
as a chemical clock, increasing their values with the evolution of 
the clumps, from intermediate, to active, to red clumps. It is not 
clea r why this trend does not include quiescent clumps (see Sec- 
tion 14.2.21 . Median values of the abundance ratios are given in 
Tabled 

we have assumed a constant value of 50 for the [ 12 C/ 13 C] 
isotopic abundance ratio for all sources, regardless of 
their different distances to the Galactic center. How- 
ever, it has been shown that there is a gradient of 
this ratio increasin g outwards from the Galactic center 
([Savage et al.ll2u0l . The [ 12 C/ 13 C] ratio has been mea- 
sured using several methods and ranges from ^20 to ^70 
([Savage et al.l [2002. and references therein). The varia- 
tion of the [ 12 C/ 13 C] ratio can give an error of a factor 
~2 in the column density estimation. 

Dust emission properties. H 2 column densities 
were calc ulated using the s moot hed 1.2 mm emission 
data from iRathborne et al.l ([2006D and the same proce- 
dure that they used. The main assumptions are that the 
dust opacity per gram of gas, Ki.2 mm , is 1.0 cm 2 g" 1 , and 
the gas-to-dust m ass ratio is 100. The value fo r Ki. 2mm 
was determined by Osscnkopf & Henning (1994) for pro- 
tostellar cores assuming dust grains with thin ice mantles 
at gas densities of 10 6 cm -3 . However, dust opacities 
and the gas-to-dust ratio are not well constrained. They 
could vary due to dust destruction in later stages of evo- 
lution (by shocks or UV radiation), or grain growth in 
dense, cold cores at very early stages. In the first sce- 
nario, the gas-to-dust ratio would be greater than 100; in 
the second scenario, the gas-to-dust ratio would be less 
than 100 and Ki. 2 mm would decrease. 

4.4. Molecules 

In this section, we discuss each commonly detected 
molecular line in depth. We describe the structure of 
the lines with hypcrfinc transitions (N 2 H + , HCN, and 
C2H), a brief review about their formation and previous 
works, and their behavior in IRDC clumps. 

4.4.1. A r 2-ff + (Diazenylium) 



Figure 17. Histograms of the number distributions of 
N2H+/HNC abundance ratios for each evolutionary sequence. The 
name of the evolutionary stage is given on the top left corner of 
each panel. The vertical dashed lines indicate the median values of 
the abundance ratios for each distribution. The diagonal patterns 
indicate the portions of the histograms that correspond to upper 
and lower limits (see bottom panel, right side). The median values 
and the K-S test suggest that the N2H+/HNC ratio may act as 
chemical clock, increasing their median values with the evolution 
of the clumps, from quiescent to red clumps. Median values of the 
abundance ratios are given in Table [8] 

The N 2 H + J = 1 — >• rotational transition exhibits a 
quadrupolc hyperfinc structure due to the spin of the two 
nitrogen nuclei. Although for N 2 H + J = 1 -4 the re are 
15 allowed hyperfme transitions (Da niel et all2006l) . only 
seven feature s can be resolved in low-mass star for ming 
regions (e.g., ICaselli et al.l [19951: ICransi et aljl2003) and 
three in high-mass star forming regions (e.g. JFoster et al.1 
2011). In the latter case, because massive star- forming 
regions show turbulent line widths (>3 km s _1 ) much 
broader than thermal widths (^0.3 km s^ 1 ), two triplets 
are blended and observed as one line. The relative inten- 
sities in the optically thin case are given by the statisti- 
cal weights 1:3:5:7, when seven components are detected, 
and 1:3:5, when three components are detected. The rel- 
ative intensities become one in the optically thick case. 

N 2 H + is a good tracer of dense, cold gas because it 
does not deplete at low temperatures and high densi- 
ties. Because the N 2 H + spectrum show no evidence of 
self-absorbed profiles or line wings indicating outflow ac- 
tivity, this molecule is also a good tracer of dense, warm 
gas. However, some clumps have broad lines, which may 
suggest they are formed by more than one central object. 
Indeed, high angular resolution observations made with 
interferometers show the clumpy nature o f IRDCs (e.g., 
IRathborne et~aII[2007L l200l IZhang et all 120091) . N 2 H+ 
is the second most frequently detected line in this survey 
and its detection rates do not depend of the star for- 
mation activity. However, the brightness temperature of 
N 2 H + is normally weak in regions without signs of star 
formation. N 2 H + column densities and abundances in- 
crease with the evolution of the clumps. Although, in 
general, an increase of column density does not neces- 
sary mean an increase of the molecular abundance, we 
see in our sample of objects that column density and 
abundance of N 2 H + tend to increase together. The in- 
crease of column density with the evolution of the clumps 
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might be explained by the acc retion of material. For in- 
stance, iChambers et alj (|2009ft showed that, on average, 
clumps with signs of star formation have smaller sizes 
and higher de nsities than clumps wit hout apparent star 
formation and Zhang & W ang! (|2011[ ) propose that em- 
bedded protostars and protostellar cores undergo simul- 
taneous mass growth during the protostellar evolution. 
The increase in the abundance might be explained by 
the rise in the temperature which releases N2 from the 
dust grains to the gas phase allowing the formation of 
more N 2 H + with the passage of time. However, with 
the increase of temperature, CO should also be released 
from the dust grains in to the gas phase and should ef- 
ficien t ly destroy N 2 H + dLee et al.|[2004l : Uorgensen et al.1 
[200l IBusquet et ail I2011D . Why exactly N 2 H+ abun- 
dance increases with evolution is not clear. We sug- 
gest two possible explanations. First, the rates involving 
the formation and destruction of N 2 , CO, N 2 H + , and 
HCO + might not be accurate, and N 2 H + might not be 
destroyed as efficiently as is currently believed. Second, 
the large beam of Mopra is including emission from cold 
gas which is sorrounding the most compact regions di- 
rectly associated with the star formation processes in 
the center of the clumps. If most of the N 2 H + emis- 
sion is coming from a cold envelope, in clumps that show 
signs of star formation, we are using in the abundance 
calculation a higher temperature than the actual value. 
Higher tempe ratures would l e ad us to infer larger N 2 H + 
abundances. IBusquet et al.l ()2011h performed a time- 
dependent chemical modeling of the massive protostel- 
lar clu ster AFGL5142 based on the model of iViti et al.l 
(2004). Their model consists of a prestellar and a pro- 
tostellar phase. From their Figure 9, it can be seen that 



once temperature increases, from the prestellar to the 
protostellar phase, the N 2 H + molecular abundance also 
increases, supporting what we see in our work. 

4.4.2. HCO+ andH vi CO + (Formylmm) 

HCO + has been widely used t o investigat e infa ll 
motions (e.g., iFuller et al.1 [20051 : iChen et all [2010h . 
and, occasionally, its spectrum shows high-velocity 
wing emission which i ndicates outflow activity (e.g., 
iCveanowski et al.l 120 111) . This makes HCO+ a good 
tracer of kinematics in star- forming regions. On the other 
hand, H 13 CO + is optically thin and normally shows a 
Gaussian profile. The combination of these two lines 
is usually used to distinguish between an asymmetric 
"blue profile" (indicative of infall motions) and two ve- 
locity components in the line of sight. If the optically 
thin H 13 CO + peaks in the self- absorption of the op- 
tically thick HCO + , it indicates a genuine "blue pro- 
file." In our sample of IRDC clumps, three intermedi- 
ate clumps (G024.33 MM2, G024.33 MM5, and G031.97 
MM5) present this characteristic profile associated with 
large-scale infall motions. As also happpens for N 2 H + , 
HCO + column densities and abundances rise with the 
evolution of the clumps. As we mentioned before for 
N 2 H + , the increase in HCO + column densities might be 
due to the physical evolution of the clumps. The rise in 
HCO + abundance can be plausibly explained by the in- 
crease in the temperature, which releases CO from grain 
mantles into the gas phase, forming more HCO + with the 
passag e of time. For example, although IBusquet et al.1 
(|201 If ) focused on the NHa/N 2 H + abundance ratio, they 
also obtained HCO + abundances from their chemical 
modeling. They showed that once temperature rises, 
from the prestellar to the protostellar phase, the HCO + 
molecular abundance rises, as well as that of N 2 H + . 



Table 8 

Median Values of Derived Clump Properties for each Evolutionary Stage 



Molecule 



Quiescent 
Clumps 



Intermediate Active 
Clumps Clumps 



Red 
Clumps 



All 
Clumps 



Optical Depths (r) 



HNC 


24 




22 


18 




16 




19 




N 2 H+ a 


0.8 




0.8 


0.5 




1.8 




0.8 




HCO+ 


30 




20 


16 




19 




21 




C 2 H b 


2.9 




2.4 


2.2 




1.4 




2.4 










Column Densities (cm 1 


) 










HNC 


1.84 x 10 


14 


2.50 X 10 14 


3.39 x 10 


14 


2.45 x 10 


14 


2.42 x 10 


14 


N 2 H+ 


1.03 x 10 


13 


1.40 X 10 13 


2.34 x 10 


13 


2.75 x 10 


1.:-! 


1.60 x 10 


13 


HCO+ 


1.40 x 10 


14 


1.51 x 10 14 


2.51 x 10 


14 


3.20 x 10 


14 


1.88 x 10 


14 


C 2 H 


2.02 x 10 


14 


2.40 x 10 14 


2.66 x 10 


14 


2.68 x 10 


14 


2.41 x 10 


14 


HC3N 


3.05 x 10 


12 


4.10 x 10 12 


5.47 x 10 


12 


7.97 x 10 


12 


4.77 x 10 


12 


HNCO 


1.33 x 10 


13 


3.36 x 10 13 


4.30 x 10 


13 


4.31 x 10 


13 


3.36 x 10 


13 


SiO 


1.36 x 10 


12 




1.39 x 10 


13 


4.62 x 10 


12 


7.72 x 10 


12 








Molecular 


Abundances 












HNC 


3.52 x 10 


-8 


4.41 X 10~ s 


3.29 X 10 


-8 


6.11 x 10 


-8 


3.73 x 10 


-8 


N 2 H+ 


1.65 x 10 


-9 


2.73 X 10" 9 


3.13 x 10 


-9 


3.73 x 10 


-9 


2.40 x 10 


-9 


HCO+ 


2.40 x 10 


-8 


2.48 X 10" 8 


4.23 x 10 


-8 


5.65 x 10 


-8 


2.51 x 10 


-8 


C 2 H 


3.65 x 10 


-8 


3.45 x 10~ 8 


4.26 x 10 


-8 


4.92 x 10 


-8 


3.72 x 10 


-8 


HC3N 


3.14 x 10~ 


-10 


2.34 x 10~ 10 


4.97 x 10" 


-10 


4.23 x 10" 


-10 


4.23 x 10- 


-10 


HNCO 


1.83 x 10 


-9 


3.21 x 10~ 9 


2.51 x 10 


-9 


4.00 x 10 


-9 


2.80 x 10 


-9 


SiO 


2.78 x 10~ 


-1(1 




5.95 x 10- 


-1(1 


4.37 x 10" 


-10 


5.79 x 10- 


-10 


Abundance Ratios 


N 2 H+/HCO+ 


0.08 




0.12 


0.10 




0.07 




0.09 




N 2 H+/HNC 


0.06 




0.06 


0.08 




0.08 




0.07 





Optical depth for the brightest transition JF\F ■ 
Optical depth for the brightest transition NJF = 



■■ 123 - 
1«2 



012. 

> 04 
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Table 9 

Column Density Comparison with Other Works. 



Different Works 


N 2 H+ 
(xlO 13 ) 
(cm -2 ) 


HCO+ 
(xlO 14 ) 
(cm -2 ) 


HNC 
(xlO 14 ) 
(cm -2 ) 


C 2 H 
(xlO 14 ) 
(cm -2 ) 


HC 3 N 
(xlO 12 ) 
(cm -2 ) 


HNCO 
(xlO 13 ) 
(cm -2 ) 


SiO 

(xlO 13 ) 
(cm -2 ) 


This Work 


Q: 1.03 
A: 2.34 


A: 2.51 
R: 3.20 


A: 3.39 


A: 2.66 


Q: 3.05 


A: 4.30 
R: 4.31 


A: 1.39 
R: 0.46 


Sakai ct al. (2008) 
Piroeov et al. (20031 
Sakai et al~T2010) 
Purcell et al"72006) 
Zinchenko et al. (2000) 
Miettinen et al. (20061 


Q: 1.15 
A: 2.78 


A: 3.30 
R: 11.7 


A: 3.60 


A: 3.25 


Q: 9.70 


A-R: 8.61 


A: 1.46 
R: 2.00 



Note. — Q — Quiescent; A — Active; R — Red. 
This table displays median values for comparison. 

4.4.3. HNC and HN 13 C (Hydrogen Isocyanide) 

HNC is the molecule with the highest detection rate in 
the sample and its presence is independent of the star for- 
mation activity. Because HNC is ubiquitous in the IRDC 
clumps, and its spectrum shows no evidence of line wings 
indicating outflows and only a few self-absorbed profiles, 
HNC seems to be a good tracer of cold and warm gas. 
The HNC emission lines show, in general, broad widths 
(as seen in the HCO + lines as well, Figure [SJ. This is 
probably because both lines are optically thick, and AV 
can be broader by opacity. On the other hand, HN 13 C 
is optically thin and normally shows a Gaussian profile. 
We note that the assumption of HNC emission being op- 
tically thin produces column densities of at least one or- 
der of magnitude lower. There is a slight increase of HNC 
column densities with the evolution of the clumps. How- 
ever, there are no differences in the HNC abundances 
for the different evolutionary states. The accretion of 
material in the star-forming clumps could explain the 
small rise in the HNC column densities. From the his- 
togram with the relative abundances of N2H 4 " and HNC 
(Figure H71) , it can be inferred that HNC may be pref- 
erentially formed in cold gas. Future chemical modeling 
is needed to clarify the behavior of HNC in high-mass 
star-forming regions. 

4.4.4. HCN (Hydrogen Cyanide) 

The HCN J = 1 — > rotational transitions show three 
hyperfme components caused by the nuclear spin of 
the nitrogen nucleus. In the optically thin case, these 
three components have relative intensities of 1:3:5. How- 
ever, due to the broad linewithds found in massive star- 
forming regions, the components are blended. In addi- 
tion, several sources exhibit spectra with extended wing 
emission and a wide variety of relative intensities other 
than those in the optically thin limit. The combination of 
these three factors makes it difficult to perform Gaussian 
fits and analyze the HCN line towards the IRDC clumps 
in this survey. For a better understandi ng of HCN hyper- 
fine tr ansitions, see the recent work of iLoughnane et al.l 
(120121) . 

4.4.5. C 2 H (Ethynyl) 

C2H N = 1 — > exhibits 6 fine and hyperfme compo- 
nents due to the presence of both an electron spin and 



a nuclear spin of the hydrogen nucleus. The spectral 
separation makes all of them observable in star-forming 
regions. However, two of them are relatively very weak 
and were not detected in this survey. Their relative inten- 
sities in th e optically thin case ar e given bv lTucker et al.l 
(fl97l and lPadovani et al.l (pool . 

C2H has been known to be a PDR t r acer (e.g., 
iFuente et al"1ll993D . However. iBeuther et al.l (|2008D find 
that this molecule also seems to trace dense gas in early 
stages of star formation. They found that the distribu- 
tion of C2H shows a hole around a hot core, and suggest 
that C2H decreases in the hot core phase. Based on the 
results of their simple chemical modeling, they suggest 
that C2H may be a suitable tracer of early stages of star 
formation. 

So far, C2H has not been systematically studied in 
massive star-forming regions or IRDCs. Although we 
cannot know the spatial distribution of C2H in the IRDC 
clumps, we can say that there is definitely an increase 
in C2H detections from quiescent to active sources (see 
Figure [7]). The ~30% detection rate of this line in qui- 
escent clumps shows that its emission is not ubiquitous 
in these kinds of sources as, for example, N2H" 1 ", HNC, 
and HCO + are. Certainly, higher angular resolution ob- 
servations and mapping are needed to clarify if the C2H 
emission comes from the external layers of clumps that 
can interact with PDR emission surrounding them, or if 
the C2H emission comes from the dense, cold gas inside 
quiescent clumps. We note that C2H lines show the best 
AV correlation with N2H + , suggesting that the emis- 
sion originates from the same region (see Figure O. C2H 
column densities and abundances show no clear trend 
of changing with evolution. There is marginal evidence 
for increasing column densities with the evolution of the 
clumps, but it is not supported by the K-S test. 

4.4.6. HC3N (Cyanoacetylene) 

HC3N is the simplest of the cyanopolyynes, molecules 
in the form of HC2 n +iN with n from 1 to 5. Its main 
progenitor, C2H2, exists on grains mantles and is re- 
leased into the gas phase during the onset of heating. 
Thus, HC3N is associated with warm, dense gas in re- 
gions with current sta r formation, such as hot cores (e.g., 
iChapman etld1l2T)0l . This may explain the high detec- 
tion rates (^50%) for this molecule in active clumps. 
The presence of HC3N and the rich molecular line spec- 
tra found in active cores (see Figure [6]), suggest that 
some active clumps may have embedded hot molecular 
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cores, as, for example. iRathborne et all (|2008l ) found in H 13 CO + and HN 13 C (see Figure IH]). HC 3 N column den- 
one IRDC clump using high angular resolution observa- sities vary with the evolution of the clumps, increasing in 
tions. The HC3N linewidths are more narrow than N2H + more evolved regions as is shown by the median values, 
linewidths, and seem to correlate better with those of For HC3N abundances, no real change is clear. 



Table 10 

Continuum Emission Parameters and Molecular Abundances 



IRDC Clump 


1 2 mm 








Molecular Abundance 












r lux 


y in22 






HNC 


C 2 H 












(mjy) 


(cm - ^ ) 


x 10 — ^ 


X 10 — 8 


xlO" 8 


xlO -8 


X 1U 


x 10 — ^ 


X 10 


G015.05 


MM1 


45 


0.86 


1.82(0.15) 


2.39(0.85) 


2.64(1.27) 


3.20(3.52) 




1.94(0.81) 




G015.05 


MM2 


25 


0.35 


2.05(0.45) 


1.69(0.60) 


7.70(5.53) 










G015.05 


MM3 


18 


0.33 


1.08(0.31) 




3.12(1.52) 










G015.05 


MM4 


14 


0.18 


5.51(0.63) 


3.09(1.07) 


8.26(3.43) 


16.20(17.41) 








G015.05 


MM5 


24 


0.32 


1.63(0.51) 


3.16(1.74) 


7.50(4.84) 










G015.31 


MM2 


22 


0.43 


2.73(2.89) 




5.05(4.12) 










G015.31 


MM3 


17 


0.32 


0.96(0.25) 


2.65(1.90) 


1.58(1.61) 










G015.31 


MM5 


22 


0.34 




3.22(1.83) 


3.98(2.79) 










G018.82 


MM2 


77 


0.82 


2.74(2.75) 


3.91(2.38) 


4.21(2.55) 










G018.82 


MM3 


30 


0.36 


16.75(16.99) 


27.86(14.46) 


18.52(8.67) 










G018.82 


MM4 


■13 


1.26 


0.85(0.14) 


6.73(5.33) 


3.98(2.25) 






0.77(0.83) 




G018.82 


MM6 


46 


0.88 


1.39(0.14) 


0.60(0.21) 


0.91(0.40) 










G019.27 


MM2 


16 


0.65 


3.36(0.22) 




6.30(2.15) 


6.40(7.08) 








G022.35 


MM1 


93 


2.19 


0.19(0.05) 


0.39(0.15) 


0.51(0.21) 










G022.35 


MM2 


32 


0.29 






6.94(3.70) 










G023.60 


MM7 


19 


0.43 


9.45(9.31) 




8.00(3.82) 










G023.60 


MM9 


32 


0.41 


1.67(0.27) 




3.58(1.72) 


3.68(5.28) 








G024.08 


MM2 


18 


0.72 


0.98(0.16) 


2.49(1.01) 


3.52(1.57) 










G024.08 


MM3 


33 


0.48 


1.94(1.92) 


3.72(1.85) 


3.98(1.77) 










G024.08 


MM4 


15 


0.63 


0.76(0.21) 


3.86(3.50) 


2.92(1.49) 










G024.33 


MM2 


11 


0.56 


5.44(0.32) 


2.43(0.63) 


10.94(2.07) 






7.58(1.95) 




G024.33 


MM3 


53 


0.70 


3.77(0.24) 


1.97(0.50) 


1.22(0.41) 


4.85(6.65) 




4.01(1.28) 


5.64(2.02) 


G024.33 


MM4 


53 


0.72 


2.19(0.19) 


1.10(0.42) 


1.10(0.40) 










G024.33 


MM5 


51 


0.67 


3.24(0.22) 


4.18(0.68) 


1.18(0.42) 


3.92(5.95) 




8.10(1.39) 




G024.33 


MM7 


26 


0.49 


1.89(0.17) 


2.49(0.47) 


2.15(0.78) 


1.24(2.46) 






2.78(1.51) 


G024.33 


MM8 


12 


0.55 




2.50(0.95) 


4.69(1.80) 










G024.33 


MM9 


13 


0.31 




7.24(2.81) 


10.96(4.29) 










G024.33 


MM11 


32 


0.55 


2.76(0.17) 


5.16(1.14) 


7.34(2.73) 










G024.60 


MM2 


78 


1.54 


0.81(0.06) 


0.98(0.32) 


0.96(0.41) 


0.75(0.26) 


2.25(0.80) 






G025.04 


MM2 


59 


0.91 


0.68(0.12) 


1.39(0.56) 


1.87(0.80) 










G025.04 


MM4 


62 


0.89 


2.74(0.12) 


1.69(0.61) 


4.82(1.13) 


1.22(2.49) 




3.50(0.73) 




G027.75 


MM2 


17 


0.32 


3.12(3.15) 














G027.94 


MM1 


11 


0.37 


4.44(4.46) 




5.28(1.83) 


5.19(8.39) 








G028.04 


MM1 


11 


0.47 


4.87(4.79) 


4.80(2.04) 


11.79(5.69) 


2.73(4.87) 


8.07(2.62) 


6.17(1.80) 




G028.08 


MM1 


29 


0.61 


1.03(1.02) 




3.14(1.66) 










G028.23 


MM1 


59 


1.23 


2.54(2.55) 


0.81(0.27) 


0.62(0.22) 










G028.28 


MM4 


13 


0.51 


2.96(0.26) 


3.55(0.73) 


3.20(1.17) 


5.19(5.99) 








G028.37 


MM1 


318 


4.02 


1.14(0.04) 


0.62(0.10) 


1.42(0.24) 


1.07(0.60) 


3.60(0.36) 


2.48(0.22) 


3.98(0.34) 


G028.37 


MM2 


118 


1.23 


2.04(0.15) 


1.97(0.36) 


5.18(1.43) 


3.70(5.25) 


3.96(1.73) 


2.94(1.22) 




G028.37 


MM4 


115 


1.51 


3.10(0.12) 


1.92(0.20) 


2.89(0.31) 


4.17(3.73) 


5.65(0.71) 


5.33(0.42) 


9.20(0.65) 


G028.37 


MM6 


87 


1.71 


1.56(0.06) 


1.59(0.30) 


1.68(0.27) 


1.06(1.38) 


4.48(0.71) 


2.51(0.30) 


5.96(0.52) 


G028.37 


MM9 


67 


1.21 


1.26(0.09) 


1.64(0.29) 


1.26(0.24) 










G028.37 


MM11 


53 


0.60 


1.74(0.34) 


1.36(0.47) 


11.93(6.30) 


4.36(6.98) 








G028.37 


MM12 


51 


0.97 


0.92(0.93) 


2.49(1.42) 


3.71(2.94) 




3.14(1.58) 






G028.37 


MM13 


20 


0.28 


5.16(5.59) 














G028.53 


MM3 


85 


1.78 


1.39(1.37) 


1.19(0.36) 


1.34(0.58) 










G028.53 


MM5 


11 


0.63 


4.63(4.83) 


2.39(1.62) 


5.48(3.21) 










pniiQ erg 

0028.53 


MM / 


11 


0.86 


1.28(0.14) 


1.06(0.31) 


4.03(1.41 ) 










G028.53 


MM8 


35 


0.54 


1.84(0.26) 


6.53(5.13) 


6.35(3.27) 


3.73(2.12) 








G028.53 


MM9 


61 


0.94 


2.00(0.16) 


2.53(1.01) 


3.27(1.07) 










G028.53 


MM10 


90 


1.28 


1.59(0.15) 


2.22(0.85) 


3.83(1.14) 


3.96(5.40) 








G028.67 


MM1 


27 


0.27 




6.07(3.02) 


10.93(5.20) 










G028.67 


MM2 


38 


0.76 


0.95(0.14) 




3.51(2.09) 










G030.14 


MM1 


10 


0.64 


2.12(2.33) 


1.83(1.05) 


2.50(1.28) 










G030.57 


MM1 


117 


1.40 


3.16(0.11) 




2.75(0.87) 


1.89(2.95) 


6.10(1.11) 


2.04(0.70) 




G030.57 


MM3 


17 


0.34 


2.92(0.55) 




18.84(10.70) 


7.12(10.94) 








G030.97 


MM1 


98 


1.08 


3.38(0.14) 


4.80(0.46) 


3.75(0.62) 




3.43(1.13) 






G031.97 


MM5 


29 


0.35 


6.77(0.39) 


3.63(0.85) 


10.27(3.34) 


8.52(11.04) 








G031.97 


MM7 


36 


0.59 


2.89(0.16) 


2.05(0.40) 


5.13(4.79) 


3.66(5.80) 




1.70(0.96) 




G031.97 


MM8 


38 


0.31 


10.84(0.62) 




10.76(3.73) 










G033.69 


MM1 


97 


1.71 


2.08(0.06) 


1.00(0.15) 


1.72(0.20) 




5.61(0.75) 


2.66(0.35) 


3.09(0.60) 


G033.69 


MM2 


75 


0.73 


4.34(0.29) 


6.17(2.51) 


6.11(1.56) 


4.84(0.67) 


13.77(2.24) 


5.87(1.66) 




G033.69 


MM3 


52 


0.48 


4.97(0.58) 


7.48(3.12) 


10.73(4.33) 










G033.69 


MM4 


65 


1.00 


3.06(0.21) 


1.77(0.65) 


3.38(0.83) 




5.59(1.97) 


3.09(1.02) 




G033.69 


MM5 


44 


0.49 


4.42(0.50) 


7.90(3.93) 


8.69(3.61) 


7.70(6.26) 
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IRDC Clump 1.2 mm N(H 2 ) Molecular Abundance 

Flux xlO 22 N 2 H+ HCO+ HNC C 2 H HC 3 N HNCO SiO 

(mjy) (cm" 2 ) xlO -9 xlO -8 xlO -8 xlO -8 xlCT 10 xl0~ 9 xlO~ 10 



G033.69 MM11 


37 


0.65 






3.88(1.46) 


6.14 


2.33) 


2.28(4.42) 




G034.43 MM1 


548 


5.85 


2.34 


0.04) 


0.82(0.09) 


2.33 


0.27) 


1.32(0.81) 


5.70(0.36) 


G034.43 MM5 


89 


1.75 


1.70 


0.05) 


1.10(0.15) 


2.99 


0.45) 




3.07(0.71) 


G034.43 MM7 


59 


0.85 


1.31 


0.12) 


1.26(0.39) 


2.69 


0.70) 


2.41(3.42) 


2.43(1.38) 


G034.43 MM8 


69 


0.64 


2.39 


0.26) 


4.23(1.89) 


3.98 


1.58) 


5.56(9.33) 




G034.77 MM1 


111 


1.11 


3.68 


3.65) 


4.08(0.71) 


2.20 


0.72) 




4.21(1.17) 


G034.77 MM3 


18 


0.25 


3.55 


0.60) 


2.03(0.71) 


11.20(4.72) 


12.12(20.78) 




G035.39 MM7 


50 


0.61 


3.92 


0.27) 


9.77(1.54) 


6.21 


1.34) 


7.13(6.73) 


5.27(2.41) 


G035.59 MM1 


22 


0.25 


11.25 


12.28) 


17.36(8.70) 


7.02 


4.36) 






G035.59 MM2 


19 


0.35 


2.03 


0.29) 


5.09(1.59) 


3.05 


1.14) 






G035.59 MM3 


30 


0.42 


2.40 


0.55) 


2.69(0.70) 


3.00 


1.13) 


4.37(5.26) 


3.32(2.99) 


G036.67 MM1 


38 


0.90 


1.10 


1.13) 


1.75(1.49) 


1.13 


0.66) 


1.75(1.74) 




G036.67 MM2 


38 


1.22 






0.54(0.24) 


0.95 


0.39) 






G038.95 MM1 


104 


3.33 


0.40 


0.02) 


1.28(0.20) 


0.51 


0.15) 


0.72(0.58) 


1.41(0.58) 


G038.95 MM2 


82 


0.82 


2.10 


0.25) 


5.66(1.29) 


2.99 


1.12) 




4.23(2.10) 


G038.95 MM3 


57 


0.88 


3.40 


3.38) 


3.88(1.00) 


2.69 


1.06) 


3.20(0.46) 




G038.95 MM4 


41 


0.40 


4.89 


4.95) 


4.54(1.69) 


5.81 


2.28) 


5.36(0.97) 




G053.ll MM1 


227 


2.43 


1.28 


0.10) 


1.95(0.20) 


0.81 


0.24) 




3.28(0.57) 


G053.ll MM2 


55 


0.63 


4.39 


4.42) 


4.47(0.90) 


2.97 


0.85) 


2.58(0.98) 




G053.ll MM4 


34 


0.34 


3.29 


0.50) 


11.95(1.79) 


4.29 


1.14) 


2.79(1.06) 


4.69(3.59) 


G053.ll MM5 


25 


0.22 


13.32 


13.61) 


5.79(2.17) 


7.07 


2.69) 


4.97(1.14) 




G053.25 MM1 


52 


0.61 


2.36 


0.24) 


2.28(0.75) 


1.98 


0.68) 


2.54(3.40) 




G053.25 MM3 


17 


0.20 


6.87 


7.23) 


7.48(3.09) 


6.77 


2.63) 


15.56(26.01) 




G053.25 MM4 


37 


0.40 


3.99 


0.37) 


6.01(1.37) 


3.04 


1.07) 






G053.25 MM5 


17 


0.20 




12.49(4.32) 


6.76 


2.55) 


10.04(22.70) 




G053.25 MM6 


34 


0.30 


13.05 


13.20) 


8.97(2.12) 


5.61 


2.10) 






G053.31 MM2 


34 


0.59 


2.09 


2.07) 


1.13(0.42) 


1.69 


0.65) 







1.47(0.25) 5.94(0.36) 



0.27(0.16) 



16(4.48) 



4.4.7. HNCO (Isocyanic Acid) 

HNCO is a high density tracer that has been found 
in regions that span a l arge range of temperatures 
(10-500 K; Uacksonet all li"98l IZinchenko etaTl [20001: 
iBisschop et al.l 120071) . The formation of HNCO is 
inefficient if on l y gas -phase reactions are considered 
(jTideswell et al.1 mm . In hot cores, the gas forma- 
tion chemistry can explain the observed HNCO abun- 
dances when reactions on the surfaces of grains are 
included. HNCO is efficiently formed on grain man- 
tles and processed in more complex molecules. How- 
ever, HNCO is not directly ejected from the dust. In- 
stead, it is formed by the dissociation of more complex 
molecules (formed by HNCO on t he grain mantles) once 
they are released to the gas phase (fTideswell et al.ll2010b 
IRodriguez-Fernandez et all 120101) . IZinchenko et all 
(2000) found evidence that HNCO and SiO may have 
a common producti on mechanism, presumably based 
on shock chemistry. IRodriguez- Fernand ez et all (j2010f ) 
found HNCO lines that exhibit the same characteristics 
as other well-known shock tracers such as CH3OH and 
some sulfur-bearing species (SO and SO2). These au- 
thors suggest that HNCO is a shock tracer which is ex- 
plained by a combination of grain surface and gas phase 
chemistry. HNCO abundances would rise in shock re- 
gions because this molecule can be directly ejected to 
the gas phase through grain sputtering and because the 
efficiency of neutral-neutral reactions increases in the gas 
phase due to higher gas temperatures. 

In our IRDC sample, HNCO emission is detected 
slightly more often in regions with signs of star forma- 



tion. The HNCO line widths seem to correlate with those 
of N2H+. HNCO profiles show no evident signatures of 
being a tracer of shocks in most of the sources, except 
in G028.37 MM1 and G028.37 MM4 (which represent 
10% of the sources with HNCO detection). In these two 
sources, the HNCO spectrum presents a blue wing which 
is also observed in SiO. HNCO column densities rise with 
the evolution of the clumps, as is shown by the median 
values. However, HNCO abundances do not show this 
trend. 

4.4.8. SiO (Silicon Monoxide) 

It is well-known that SiO is a powerful tracer of 
molecular gas asso ciated with shocks (Schilkc et al.lll997l : 
iCaselli et allll997t) . Its abundance is highly enhanced in 
molecular outflows, with respect to the ambient abun- 
dances. In_som£_cases_these enhancements can be up 
to 10 6 (jMartin-Pintado et al.lfl992l ). The increase of the 
SiO abundance is due to the silicate grain destruction, 
giving rise to the injection into gas of Si atoms and/or 
Si-bearing species, an d the subsequent hig h temperature 
gas-p hase chemistry (Schilkc ct al. 1997; Cascll i et all 
1997). The previous SiO production mechanism ex- 
plains spectral lines that show extended wing emission 
with broad line widths caused by the interaction of out- 
flows and the surrounding medium. However, an SiO 
componen t coming from quiescent g as is not well un- 
derstood. Uimenez-Serra et al.l ((2010) detected extended 
narrow SiO emission not associated with signs of star 
formation in an IRDC. They suggest that this SiO com- 
ponent could be produced by the following processes: 
remnants of large-scale shocks caused by the formation 
process of the IRDC, decelerated shocked gas associ- 
ated with large scale outflows from neighboring massive 
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protostars, recently processed material associated with 
the youngest jets/outflows, and/or an undetected and 
widespread lower mass protostar population. 

In our IRDC sample, SiO emission is mostly detected 
in clumps with signs of star formation. The only excep- 
tion is the quiescent clump G024.33 MM7 which shows 
a narrow line profile (1.4 km s _1 ) and the lowest SiO 
abundance in the sample. The average SiO line width in 
the sample is ~5 km s _1 . Half of the sources with SiO 
detection show clear extended wing emission indicating 
the presence of outflow activity. 

5. CONCLUSIONS 

We have carried out a multi-line survey at 3 mm to- 
ward 37 IRDCs, containing 159 clumps, in order to inves- 
tigate the behavior of the different molecular tracers and 
search for chemical variations through an evolutionary 
sequence based on Spitzer IRAC and MIPS emission. We 
observed N 2 H+, HNC, HN 13 C, HCO+, H 13 CO+, HCN, 
C 2 H, HC 3 N, HNCO, and SiO lines with the Mopra 22 m 
telescope located in Australia. After eliminating clumps 
that are not located in IRDCs and pairs of clumps that 
are placed within one Mopra beam, we base our study 
on 92 sources. 

HNC and N2H 4 " lines are detected in almost every 
IRDC clumps at every evolutionary stage, indicating that 
their presence does not depend on the star formation ac- 
tivity. On the other hand, HC 3 N, HNCO, and SiO lines 
are predominantly detected in later stages of evolution, 
as expected from their formation paths. 

The line widths of N2H 4 " slightly increase with the evo- 
lution of the clumps, which is likely produced by the rise 
of turbulence due to the enhancement of the star forma- 
tion activity at later evolutionary stages. The increase 
is modest because, due to the large Mopra beam, we are 
also tracing the bulk motions of the gas, instead of just 
the densest regions associated with star formation. 

Optical depth calculations show that the N2H 4 " line is 
mostly optically thin (median of 0.8) and the C2H line is 
moderately optically thick (median of 2.4). HCO + and 
HNC lines are optically thick (medians of 21 and 19, re- 
spectively), while their isotopologues are optically thin 
(median of 0.4 for both). N2H + opacities show no vari- 
ations with the evolution of the clumps, whereas C2H, 
HCO + and HNC show a slight decrease with the rise of 
star formation activity. 

In general, column densities of the different molecules 
change for the differen t evolutionary stages defined by 
iChambers et all ([2009) and increase with the evolution 
of the clumps, with the exception of C2H. However, this 
is not generally true for molecular abundances (i.e., af- 
ter dividing by the total H2 column density inferred 
from 1.2 mm continuum emission). Only the increases 
of N2H + and HCO + abundances are statistically signif- 
icant and reflect c hemical evolution. Th is is consistent 
with the results of iBusauet et al.l (|2011l ). who included 
both molecules in their chemical modeling of a massive 
star-forming region. Although it is expected a rise of the 
HCO + abundance with the evolution of the clumps, it is 
not clear why N2H + also follows this trend. 

The N2H + /HCO + abundance ratio acts as a chemi- 
cal clock, increasing its value from intermediate to ac- 
tive and red clumps. This observed trend is consistent 
with the theoretical predictions. The chemical models 



suggest that when clumps warm up, they release CO 
from grain mantles to the gas phase. This rise of CO 
increases the amount of HCO + , because CO is its main 
supplier, with respect to N2H + , because CO is its main 
destroyer. However, the observed trend does not extend 
to quiescent clumps. This could be due to observational 
limitations or because the theoretical predictions are in- 
correct at very early stages of evolution. We also find 
that the N2H+/HNC abundance ratio increases with the 
evolution of the clumps, from quiescent to red clumps. 
It is not clear why this ratio behaves in this way, but it 
suggests that HNC may be preferentially formed in cold 
gas. 
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